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for  the  analysis  which  calculates  the  mean  overall  heat 
transfer  coefficient  and  the  overall  pressure  drop  for  many 
flow  arrangements,  taking  into  account  the  variation  of  the 
heat  transfer  coefficients  and  the  pressure  drop  with  tempera¬ 
ture  and/or  length  of  flow  path. 

The  code  is  not  limited  to  surfaces  found  in  the  litera¬ 
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NOMENCLATURE 


English  Letter  Symbols 

2 

A  -  total  heat  transfer  area,  in 

c  -  specific  heat,  BTU/lbm-°F 

P 

C  -  heat  capacity  rate,  BTU/hr-°F  =  m  cp 
D  -  diameter,  in. 

f  -  friction  factor 

F  -  LMTD  correction  factor 

2 

gc  -  acceleration  of  gravity,  32.2  ft/sec 

2 

H  -  corrected  heat  transfer  coefficient,  BTU/hr-ft 
h  -  bank  height,  in. 

J  -  Colburn  factor 

k  -  thermal  conductivity,  BTU/hr-ft- °F 

L  -  length,  in. 

l  -  fin  height,  in. 

m  -  mass  flow  rate,  lbm/hr 

m  -  v'2H  /k-t  ,  ft-1 

o  f 

n  -  number  of  _  (used  with  appropriate  subscript) 

N  -  number  of  tubes 

p  -  pressure,  psi 

P  -  pitch,  in. 

Pr  -  Prandtl  number 

0  -  heat  transfer  rate,  Btu/hr 

r  -  radius,  in. 

2 

R  -  heat  transfer  resistance,  hr-ft  -°F/BTU 
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Re  -  Reynolds  number 

s  -  distance  between  adjacent  fins,  in. 

S  -  fin  spacing  center-to-center ,  in. 
t  -  fin  thickness,  in. 

T  -  true  temperature,  deg. 

U  -  overall  heat  transfer  coefficient,  BTU/hr- 
w  -  bank  width,  in. 

Greek  Letter  Symbols 

AT  -  temperature  difference,  deg. 

n  -  surface  efficiency 

<p  -  fin  efficiency 

ji  -  viscosity,  lbm/ft-hr 

o  -  density,  lbm/ft3 

4  -  temperature  correction 


Subscripts 

a 

- 

air 

c 

- 

cold 

f 

- 

fin 

ff 

- 

free  face 

h 

- 

hot 

i 

- 

inside 

j 

- 

reference  number 

i 

- 

limiting 

L 

- 

longitudinal 

m 

- 

mean 

o 

— 

outside 
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p  -  pass 

r  -  rows 

t  -  transverse 

T  -  total 

w  -  water 

x  -  cross-sectional 

1  -  entering 

2  -  leaving 

“  -  ambient 

Superscripts 

uncorrected 
*  -  estimated 

i  -  initial 


(flow) 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  cooling  of  fluids  by  passing  ambient  air  over  extended 
tube  surfaces  is  a  relatively  recent  development  in  heat 
exchangers.  Its  application  has  come  about  cautiously,  due 
to  the  usual  reluctance  to  change  from  well  established  and 
well  documented  methods,  i.e.  the  shell-and-tube  heat  exchanger 

However,  concern  for  the  environment  and  economic  pres¬ 
sures  have  necessitated  the  use  of  air  as  a  coolant.  Smith 
[1]  has  listed  some  typical  advantages  of  direct  cooling  with 
air  as  compared  to  cooling  with  water  in  a  shell-and-tube 
exchanger: 

a.  Eliminates  the  problem  of  temperature  rise  in,  and 
pollution  of ,  water  resources . 

b.  Enables  plant  location  to  be  independent  of  a  water 
supply . 

c.  Eliminates  the  necessity  of  much  coolant  piping. 

d.  Reduces  heat  exchanger  maintenance  costs  by  eliminating 
the  need  of  descaling  water-side  surfaces.  The 
mechanical  drives  will  operate  in  a  noncorrosive 
atmosphere . 

e.  Eliminates  water  treatment. 

f.  Limiting  coolant  temperatures  is  unnecessary. 

g.  Enables  installation  of  exchangers  at  elevations  above 
other  operating  equipment  at  no  penalty,  thus  reducing 
ground  area  requirements. 

With  air  colling  becoming  more  and  more  competitive  with 


water;  even  when  water  supplies  are  plentiful,  an  automated 
air-cooled  heat  exchanger  design  package  (which  could  be  used 


for  trade-off  studies,  first  cut  analysis,  and  conceptual 
design),  would  be  of  great  use. 

The  design  of  an  air-cooled  cross  flow  heat  exchanger  is 
a  complex  task  requiring  the  examination  and  optimization  of 
a  wide  variety  of  heat  transfer  surfaces.  Studies  have  shown 
that  a  poor  choice  of  either  the  heat  transfer  surfaces  or 
design  parameters  can  more  than  double  the  costs  charageable 
to  a  heat  exchanger  [2]  . 

For  the  optimized  design  of  heat  exchangers  with  the 
computer;  reliable,  but  fast,  calculation  methods  for  the 
mean  overall  heat  transfer  coefficient  and  the  overall 
pressure  drops  are  needed  for  the  following  reasons; 

a.  Conventional  simple  methods  using  mean  values  of 
temperatures  as  reference  temperatures  can  lead  to 
undesirable  errors  [3] . 

b.  Numerical  stepwise  integrations  are  prohibitively 
time  consuming . 

B.  REVIEW 

A  number  of  heat  exchanger  design  methods  have  been  proposed 
to  determine  the  optimum  heat  exchanger  design.  Bergles,  et 
al.  [4] ,  performed  an  evaluation  of  different  objective  func¬ 
tions  for  compact  heat  exchangers  with  different  heat  trans¬ 
fer  surfaces,  but  the  same  specifications.  The  method  did  not 
include  any  actual  optimization  techniques,  but  results  did 
show  that  a  great  improvement  in  heat  exchanger  performance 
can  be  made  by  proper  selection  of  design  parameters. 
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The  method  of  Fax  and  Mills  [5] ,  used  Lagrange  multi¬ 
pliers  to  optimize  a  heat  exchanger  design  under  specified 
constraints.  This  technique  required  that  the  objective 
function  and  constraints  be  expressed  explicitly  and  be  differ¬ 
entiable  throughout  the  range  of  interest.  The  total  num¬ 
ber  of  constraints  had  to  be  less  than  the  total  number  of 
variables,  and  all  constraints  had  to  be  equality  constraints. 
The  method  was  obviously  restricted  to  a  very  limited  number 
of  problems. 

For  a  unit  section  of  an  air-cooled  exchanger  of  standard 
length  and  width  equipped  with  specified  fans,  both  Schoonman 
[6] ,  and  Joyce  [7] ,  used  factorial  searches  of  fin  spacing, 
number  of  rows  and  air  rate  to  maximize  the  ratio  of  heat 
transfer  to  cost.  Because  the  number  of  sections  must  be 
rounded  to  an  integer,  this  method  only  appears  suitable  for 
large  exchangers. 

Nakayama  [3] ,  used  a  similar  approach  requiring  the  plotting 
of  heat  transfer  coefficients  and  pressure  drops  which  pre¬ 
sumably  could  be  programmed  into  a  factorial  search. 

Kern  [8]  derived  an  analytic  expression  for  the  annual 
cost  as  a  function  of  air  rate  and  number  of  rows.  (This  was 
done  by  using  ?  constant  tubeside  heat  transfer  coefficient 
and  an  arithmetic  mean  temperature  difference.)  The  optimum 
was  found  by  setting  partial  derivatives  equal  to  zero. 
Alternatively,  if  the  airside  is  assumed  to  control,  the 
optimum  allocation  of  the  components  of  total  cost  might  be 
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found  through  geometric  programming,  as  illustrated  by  Auriel 
and  Wilde  [9] .  Oshwald  and  Kochenberger  [10]  ,  also  presented 
a  geometric  programming  method  for  heat  exchanger  optimization 
and  used  it  to  select  heat  exchanger  fluids  considering  power 
requirements,  cost,  tube  diameters,  velocities,  temperature, 
and  other  physical  properties. 

After  discussing  factorial,  univariate,  and  random  search 
methods  for  an  optimum  design  of  a  shell-and-tube  exchanger, 
Briggs  and  Evans  [11,  121,  discuss  a  "logical  search  method"; 
or  what  Peters  and  Nicole  [13],  call  "heuristics".  With  this 
method,  an  engineer  makes  use  of  selecting  design  variables 
close  to  the  optimum  to  obtain  an  optimum  design.  The  heuris¬ 
tic  method  is  less  scientific,  and  is  useful  only  when  computer 
time  and  storage  are  at  a  premium.  Due  to  the  large  number 
of  discrete  variables  encountered  in  air-cooled  heat  exchanger 
design,  Peters  and  Nicole  [13]  ,  chose  to  base  their  cost- 
optimizing  design  programs  on  heuristic  algorithms  (starting 
close  to  optimum) ,  specific  to  the  equipment  under  considera¬ 
tion. 

Mott,  et  al.  [2],  discuss  a  computerized  procedure  for 
designing  a  minimum  cost  heat  exchanger.  The  method  minimizes 
a  cost  index  expressed  as  a  function  of  fluid  pumping  power. 

The  algorithm  imposed  no  constraints. 

To  this  point,  no  mention  has  been  made  of  applying  the 
concepts  and  techniques  of  nonlinear  programming  to  optimizing 
the  design  of  heat  exchangers.  However,  Palen,  et  al .  [14], 
in  1974  proposed  using  the  Complex  Method  [15]  ,  for  the  heat 
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exchanger  optimization  problem.  They  found  a  minimum  cost 
shell-and-tube  exchanger  by  varying  six  geometrical  parameters. 
The  Complex  Method  requires  several  feasible  starting  designs 
before  optimization  can  be  performed. 

Johnson,  et  al .  [16],  coupled  an  existing  shell-and-tube 
condenser  design  code  with  a  constrained  function  minimization 
code  to  produce  an  automated  marine  condenser  design  program 
of  vastly  different  complexity. 

The  most  complete  work  to  date  has  been  accomplished  by 
Afimiwala  [17].  He  has  applied  various  nonlinear  programming 
methods  of  optimization  to  the  heat  exchanger  design  problem; 
including  an  experimental  interactive  graphical  approach  and 
exterior  penalty  function  techniques.  The  gradient  based 
search  methods  of  Davidson  -  Fletcher  -  Powell  and  conjugate 
gradient  were  used  for  the  resulting  unconstrained  minimiza¬ 
tions.  The  exterior  penalty  method  is  extremely  useful,  since 
an  initial  solution  satisfying  the  constraints  is  not  required. 
The  gradient  based  search  methods  are  efficient  when  considering 
computer  time. 

Finally,  Fontein  and  Wassink  [18]  utilized  the  complex 
method  of  Nelder  and  Mead  [19],  and  a  steepest  descent  method 
[20],  for  optimizing  a  shell-and-tube  exchanger. 

It  can  be  seen  that  alt-hough  there  are  many  methods  that 
have  been  presented  for  heat  exchanger  optimization,  each  of 
the  methods  has  its  own  limitations;  none  is  completely  general. 
Of  all  the  design  procedures  cited  above  (those  of  which  are 
applicable  to  cross  flow  air-cooled  heat  exchangers) ,  all  are 
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limited  to  the  120  individual  surfaces  found  in  the  open 
literature  [21]  for  the  calculation  of  the  air-side  heat 
transfer  coefficient  and  friction  factor.  Therefore,  the 
designer  is  faced  with  choosing  an  optimum  surface  from  a 
number  of  individual  optimal  designs  calculated  from  one  of 
the  above  methods.  In  addition,  the  above  methods  treat  the 
overall  heat  transfer  coefficient  as  a  constant,  or  they  become 
involved  with  time-consuming  numerical  stepwise  integrations 
in  an  attempt  to  account  for  the  varying  heat  transfer  coeffi¬ 
cients  . 

This  paper  tries  to  bridge  this  gap  by  presenting  an 
optimization  routine  that:  selects  an  optimal  surface,  takes 
into  account  the  varying  heat  transfer  coefficients  and  fric¬ 
tion  factors  across  the  exchanger,  performs  each  analysis  in 
an  iterative-free  manner,  and  may  start  with  an  infeasible 
design. 

C .  METHODOLOGY 

With  the  Control  Program  for  Engineering  Synthesis  and 
Constrained  Function  Minimization  (COPES/CONMIN)  optimizing 
scheme,  a  nonlinear  optimization  program  is  available  that  is 
capable  of  optimizing  a  wide  class  of  engineering  problems 
[22,23].  Therefore,  for  the  heat  exchanger  design  problem, 
it  was  necessary  to  develop  a  subroutine,  which  given  a  start¬ 
ing  design,  would  analyze  an  air-cooled  heat  exchanger,  and 
which  would  be  suitable  for  coupling  with  the  optimizer, 
COPES/CONMIN. 


Figure  1  illustrates  the  procedure  by  which  the  heat 
exchanger  was  analyzed.  Initial  input  consisted  of  a  com¬ 
plete  listing  of  all  design  parameters,  whether  known  or 
estimated.  Those  that  were  estimates,  i.e.  unknowns,  were 
later  passed  to  the  optimizer  as  design  variables. 

The  analysis  scheme  and  the  optimizer  will  be  discussed 
in  much  greater  detail  in  the  following  chapters. 

D.  OBJECTIVES 

The  objectives  of  this  thesis  are  two-fold. 

The  first  objective  was  to  develop  a  computer  code, 
hereafter  referred  to  as  ANALIZ,  which  would  analyze  an 
air-cooled  heat  exchanger  given  any  initial  design.  The 
analysis  scheme  was  to:  be  iterative-free,  take  into  account 
the  variation  of  the  heat  transfer  coefficients  and  the  pres¬ 
sure  drop  with  temperature  and/or  length  of  flow  path,  and 
finally,  be  written  in  such  a  manner  that  it  could  be  coupled 
with  an  existing  optimizer  COPES/CONMIN. 

The  second  objective  was  to  actually  couple  ANALIZ  with 
a  numerical  optimization  program.  This  would  produce  a 
detailed  design  program  which  would  have  the  capability 
to  determine  an  optimum  surface,  while  optimizing  the  objec¬ 
tive  function  such  as  size,  weight,  cost,  etc. 


II.  NUMERICAL  OPTIMIZATION 


A.  BACKGROUND 

Almost  all  design  problems  require  either  the  maximiza¬ 
tion  or  the  minimization  of  some  parameter  or  function. 

This  parameter  shall  be  called  the  design's  objective  function 
[24] .  For  example,  the  problem  may  call  for  a  heat  exchanger 
with  a  minimum  volume.  The  expression  for  volume  would  be 
the  design  objective  function.  For  the  design  to  be  accepta¬ 
ble,  it  must  satisfy  certain  design  constraints.  For  exam¬ 
ple,  an  air  heater  must  be  designed  so  that  it  will  fit  into 
a  given  space.  Therefore,  the  engineer  must  set  design  con¬ 
straints  on  the  maximum  size  of  the  exchanger. 

If  the  objective  function  could  be  easily  formulated 
analytically,  the  maxima  or  minima  could  be  found  by  using 
the  methods  of  differential  calculus.  However,  the  limita¬ 
tions  of  this  method  are  obvious. 

Another  numerical  method  that  would  be  satisfactory  for 
small  scale  problems  would  be  an  iterative  solution  technique. 
A  computer  program  could  be  written  containing  a  series  of 
nested  iteration  loops  that  would  vary  the  design  parameters 
and  solve  the  problem  for  a  variety  of  values  for  each  of  the 
parameters.  For  other  than  small,  easily  formulated  problems, 
the  cost  in  central  processor  (CPU) ,  time  would  be  prohibitive. 

Over  the  last  twenty  years,  many  numerical  optimization 
techniques  have  been  developed  specifically  for  computer 


utilization.  These  techniques  usually  do  not  require  a 
specific  algebraic  equation,  but  rather  any  computer  algorithm 
to  which  design  variables  can  be  input  and  from  which  the 
objective  function  and  design  constraint  values  can  be 
determined  is  acceptable.  For  this  reason,  nonlinear  pro¬ 
gramming  methods  were  chosen  for  the  air-cooled  heat  exchanger 
design.  Some  of  these  techniques  were  summarized  by  Shah  et 
al .  [25]  in  figure  2. 

1.  One-Dimensional  Search  Methods.  Two  of  the  most  common 
of  these  search  methods  are  the  golden  section  [26] 

and  quadratic  interpolation  [27]  .  The  former  iso¬ 
lates  the  minimum  in  regions  of  successively  decreasing 
size,  the  latter  performs  a  series  of  iterations 
approximating  the  objective  function  as  a  quadratic. 

2.  Multidimensional  Unconstrained  Search  Methods.  These 
unconstrained  searches  can  be  performed  by  a  sequence 
of  one-dimensional  minimizations  in  the  proper  direc¬ 
tions  . 

3 .  Multidimensional  Constrained  Search  Methods.  A  common 
method  here  for  enforcing  the  constraints  in  an  optimi¬ 
zation  scheme  are  based  on  the  sequential  penalty  func¬ 
tion  method.  These  techniques  convert  the  constrained 
optimization  problem  into  a  sequence  of  unconstrained 
problems.  This  is  accomplished  by  applying  either  an 
exterior  or  interior  penalty  to  the  objective  function. 
The  Complex  Method  [15],  locates  the  optimum  based  on 
an  intuitive  approach  in  a  n-dimensional  space  defined 
by  the  independent  design  variables.  The  method  of 
feasible  directions  is  used  primarily  for  inequality 
constraints . 

An  optimization  program  based  upon  the  Augmented  Lagrangian 
Multiplier  Method  and  the  method  of  feasible  directions  was 
chosen  for  this  research  project. 


B.  CONSTRAINED  FUNCTION  MINIMIZATION  (CONMIN) 

Vanderplaats  [22],  developed  an  optimization  program, 
CONMIN,  based  on  the  method  of  feasible  directions  which  is 
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capable  of  optimizing  a  wide  variety  of  engineering  problems. 
CONMIN  is  a  FORTRAN  program,  in  subprogram  form,  that  opti¬ 
mizes  a  function  subject  to  a  set  of  inequality  constraints. 

The  following  definitions  will  be  useful  in  the  following 
discussion : 

1.  Design  Variables  -  those  parameters  which  the  optimiza¬ 
tion  program  can  change  in  order  to  improve  the  design. 

2.  Design  Constraints  -  those  parameters  which  must  not 
exceed  given  bounds  for  the  design  to  be  acceptable. 

3.  Objective  Function  -  the  parameter  which  is  going  to 
be  minimized  or  maximized. 

The  general  nonlinear  inequality  constrained  optimization 
problem  can  be  written  mathematically  as  follows  [28]: 

(1) 

0  j  »  1 , NCON  (2) 

1  k  =  1 , NDV  (3) 

(4) 


The  vector  X  is  the  vector  of  design  variables,  with  NDV 

equal  to  the  number  of  design  variables.  The  objective  function. 


Minimize  F(X) 
Subject  to: 

gj(X)  < 


x:  <  x. 

i—i 


where 


X  = 


NDV 
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F (X) ,  given  by  eq.  (1),  as  well  as  the  constraint  functions 
given  by  eq.  (2) ,  may  be  linear  or  nonlinear  functions  of 
the  design  variables.  They  shall  also  be  explicit  or  implicit 
functions  of  X,  but  must  have  continuous  first  derivatives. 

•  •  £  U 

NCON  is  the  number  of  constraints.  X^  and  X^  are  the  lower 
and  upper  bounds  or  side  constraints  placed  on  the  design 
variables.  Side  constraints  could  be  included  in  eq.  (2), 
but  are  treated  separately  for  efficiency.  Equality  con¬ 
straints  are  not  dealt  with  by  CONMIN,  but  will  be  treated 
separately  by  a  multiplier  method. 

CONMIN  requires  that  an  initial  vector  of  design  variables, 

X,  which  may  or  may  not  yield  a  feasible  design,  be  specified. 
The  design  process  continues  iteratively  as: 

X<3+1  =  X*3  +  a*  S*3  (5) 

where  S*3  is  a  vector  search  direction,  ct*  is  a  scalar  quantity 
which  determines  the  amount  of  change  in  X  and  q  is  the  itera¬ 
tion  number.  At  iteration  q  a  direction  S*3  must  be  found 
which  will  reduce  the  objective  (usable  direction  for  mini¬ 
mization)  ,  without  violating  any  constraints  (feasible  direction', 
see  Figure  3  [30]  .  Once  S*3  is  determined,  eq.  (5)  becomes  a 
one-dimensional  search  problem  in  whcih  a*  must  be  found 
such  that*  F(X)  is  at  a  minimum  (see  Figure  4  [29]),  a  new 
constraint  is  encountered,  or  a  currently  active  constraint 
(g^ (X)  =  0)  is  encountered  again. 

The  design  problem  at  iteration  q+1  becomes  one  of 
finding  a  usable-feasible  direction,  sf3,  and  a  move  parameter 


a*.  This  process  is  illustrated  geometrically  by  Johnson 
[30] .  Consider  a  condenser  problem  with  just  two  design 
variables,  and  X2,  where 

X^  =  condenser  tube  outside  diameter, 

X2  =  tube  pitch  to  diameter  ratio. 

Let  the  objective  function  be  condenser  volume,  VOL  (X) . 

Assume  that  the  tube  bundle  diameter  must  be  greater  than  a 
given  value,  BDmin/  and  that  the  cooling  water  pumping  power 
must  be  less  than  a  given  value  HP  .  Figure  5  illustrates 
the  design  problem  geometrically. 

It  should  be  reiterated  here,  that  while  Johnson's  exam¬ 
ple  starts  with  a  feasible  initial  design,  @  ,  this  is  not 

a  requirement  and  CONMIN  is  capable  of  optimizing  given  an 
infeasible  initial  design.  This  is  obviously  a  great  benefit. 

The  optimization  process  begins  by  calculating  the  gradient 
of  the  objective  function  by  using  finite  difference.  Each 
design  variable  is  perturbed  by  .01  in  a  single  forward  step. 

The  gradient  of  the  objective  function,  VF ,  shown  in 
Figure  6  is  simply  the  vector  of  the  first  partial  deriva¬ 
tives  with  respect  to  the  design  variables;  that  is: 


► 


Therefore,  because  no  constraints  are  active  or  violated  at 
(X)  .  VF  defines  the  direction  of  steepest  ascent.  Because 
it  is  desired  to  minimize  F,  the  greatest  improvement  can  be 
made  by  moving  in  the  negative  gradient  direction  so  that 


S  =  -  VF  =  -  7V0L 


(7) 


With  the  value  of'  S  now  determined,  a  search  is  per¬ 
formed  from  until  the  minimum  F  is  found  at  on  Figure 

6 .  This  is  accomplished  by  taking  several  values  of  in 

eg.  (5)  and  interpolating  for  the  a*  which  will  give  the  mini¬ 
mum  value  of  F . 

The  second  design  iteration  is  begun  at  by  again 
perturbing  X  to  find  VF .  Instead  of  moving  in  the  direction 
of  steepest  descent,  a  conjugate  direction,  developed  by 
Fletcher  and  Reeves  (31]  ,  is  chosen  with  this  method,  S  is 
calculated  as  follows: 


S*3  = 


-  VF (X)q  + 


I  VF(X)q  1  2 
|VF(X)q  r|2 


s^1  , 


(8) 


see  Figure  7. 

The  Fletcher-Reeves  method  is  used  in  order  to  speed 
convergence.  With  the  new  conjugate  direction,  a  search  is 
performed  in  this  direction  until  a  constraint  is  encountered. 
This  occurs  at  of  Figure  3  on  the  pumping  power  constraint. 

At  ©  ,  with  the  HP  constraint  active,  not  only  is  VF 
found,  but  the  gradient  of  the  active  constraint  is  also 
computed,  again  using  finite  difference.  The  requirements  on 
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the  new  search  direction  are  now  twofold;  it  must  reduce  the 
objective  function  and,  at  the  same  time,  not  violate  the 
active  constraint.  This  is  solved  by  using  the  method  of 
feasible  directions  developed  by  Zoutendijk  [32] ,  and  imple¬ 
mented  by  Vanderplaats  and  Moses  [33] . 

The  problem  of  finding  the  new  S  can  be  stated  as  [34] : 

Maximize  6  subject  to  the  constraints 


7F(X)  *  S  +  3  <_  0 

(9) 

7g j (X) *S  +  3.3  1  0 

j  =  1 , NAC 

(10) 

S  •  S  <  1 

(ID 

where  Vg^ (X)  =  -  VHP  and  NAC  is  the  number  of  active  constraints 
(in  this  case  NAC  =  1)  . 

If  equation  (9)  is  satisfied  and  6  is  positive,  the  search 
direction  will  reduce  the  objective  function  and  is  defined 
as  a  usable.  If  equation  (10)  is  satisfied  and  3  is  positive, 

S  is  a  feasible  direction,  because  in  this  direction,  no  con¬ 
straints  will  be  violated  if  only  a  small  move  is  taken.  See 
Figure  3.  0 ^  is  defined  as  the  push-off  factor  for  the  active 

constraint  and  causes  the  design  to  move  away  from  the  con¬ 
straint.  Gj  must  be  greater  than  or  equal  to  zero  in  order 
to  maintain  a  feasible  design.  If  the  maximum  value  of  3  from 
equations  (9)  through  (11)  is  zero,  then  there  is  no  direction 
that  will  both  reduce  the  objective  function  and  also  be  feasi¬ 
ble.  Therefore,  the  current  design  is,  at  least,  a  local 
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minimum.  In  Johnson's  example,  a  usable-feasible  direction 

exists  and  a  cne-dimensional  search  leads  to  in  Figure  8 

where  the  minimum  bundle  diameter,  BD  .  ,  constraint  is  met. 

'  min' 

From  it  should  be  noted  here  that  CONMIN  had  informa¬ 
tion  regarding  the  linearity  of  the  constraint  and, 

therefore,  in  (10)  ,  has  set  9^  =  0  to  allow  S  to  follow  the 
constraint  as  shown  in  Figure  9.  The  one-dimensional  search 
along  this  constraint  is  carried  out  until  no  further  design 
improvement  is  realized.  This  occurs  at  . 

This  discussion  of  CONMIN  would  not  be  complete  without 
citing  the  program's  limitations.  NDV  directly  affects  the 
computational  time  required  to  reach  the  optimum.  Since  the 
calculation  of  gradients  required  for  each  design  variable  at 
the  beginning  of  each  design  iteration  is  found  by  a  finite 
difference  step,  which  requires  a  complete  pass  through  the 
analysis  portion  of  the  program,  there  is  a  subsequent  in¬ 
crease  in  CPU  time  as  NDV  increases.  Also,  due  to  the  inter¬ 
action  between  design  variables,  as  NDV  gets  larger,  convergence 
slows  during  the  optimization  process.  Vanderplaats  [24], 
recommends  that  for  most  problems  of  general  interest  a  prac¬ 
tical  limit  of  NDV  =  20  be  imposed.  NCON  does  not  present 
the  same  problem  because  gradient  information  is  calculated 
simultaneously  with  VF  and  then  only  if  the  constraint  is 
active  or  violated. 

CONMIN  offers  no  guarantees  that  a  global  minimum  has  been 
reached.  Therefore,  to  lend  some  assurance,  the  design  is 
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started  with  several  different  initial  vectors  until  the 
same  optimal  design  is  reached. 

Although  CONMIN  performs  very  well  with  inequality  con¬ 
straints,  equality  constraints  such  as: 

hK(X)  =  0 

cannot  be  dealt  with  directly,  but  must  be  treated  separately, 
using  a  different  method,  which  will  be  discussed  in  the 
following  sections. 

C.  CONTROL  PROGRAM  FOR  ENGINEERING  SYNTHESIS  (COPES) 

Recall  that  CONMIN  was  written  in  subroutine  form, 
Vanderplaats  [23]  ,  has  developed  a  main  program  which  greatly 
enhances  the  use  of  CONMIN. 

For  this  main  program,  COPES,  the  user  must  supply  an 
analysis  subroutine  titled  ANALIZ.  The  subroutines  CONMIN 
and  ANALIZ  are  then  used  by  COPES  to  optimize  the  objective 
function  subject  only  to  the  inequality  constraints. 

ANALIZ  must  be  organized  into  three  segments:  input, 
analysis  and  output.  Based  on  the  value  of  a  counter,  ICALC, 
ANALIZ  performs  the  proper  function  in  sequence. 

The  COPES  program  currently  provides  four  specific  capa¬ 
bilities: 

1.  Single  Analysis  -  one  cycle  through  the  program,  as 
if  ANALIZ  was  executing  alone. 

2.  Optimization  -  minimization  or  maximization  of  the 
objective  function  with  constraints  and  side- 
constraints  imposed. 


3.  Sensitivity  Analysis  -  used  to  explore  the  effect  of 
changing  one  or  more  design  variables  on  one  or  more 
functions . 

4.  Two  Variable  Function  Space  -  provides  tables  of  data 
of  all  specified  combination  of  two  design  variables. 

However,  a  recent  addition  to  COPES  (still  in  the  develop¬ 
mental  stages) ,  has  put  the  Augmented  Lagrangian  Multiplier 
Method  (ALMM)  at  the  disposal  of  the  programmer.  Because  of 
its  good  rate  of  convergence  and  its  theoretical  properties, 
the  ALMM  is  preferred  for  equality  constrained  problems  [35] . 

Therefore,  COPES  will  call  the  various  subroutines  in  order 
to  optimize  the  objective  function;  subroutine  ANALIZ  for 
necessary  analysis  information,  CONMIN  for  an  optimum  based 
on  inequality  constraints,  and  a  subroutine  (yet  to  be  named), 
utilizing  ALMM  for  an  optimum  satisfying  the  equality  con¬ 
straints  . 

For  detailed  explanation  of  the  ALMM,  its  background  and 
mathematical  derivation,  consult  reference  [35];  but  for  now 
consider  the  equality  constrained  problem: 

Min  f (X)  (10a) 

Subject  to  hR(X)  =  0  k  =  1,NEC0N  (10b) 

where  NECON  is  the  number  of  equality  constraints. 

Define  the  "modified"  Lagrangian  function  as: 

_  m  __ 

L(X,\)  =  f  (X)  +  l  X.  h.  (X) 

i=l  1  * 

where  A ^  is  the  Lagrangian  multiplier.  The  problem  can  now 
be  stated  as: 


■i 


f- 
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Min  L (X , A ) 


(10c) 


Subject  to  h^(X)  =  0  (10d) 

Then,  according  to  Lagrange,  if  a  I  can  be  found  for  which 
X  solves  the  problem  stated  above,  then  X  is  also  the  solution 
to  the  original  problem,  eq.  (10a)  and  (10b) . 

The  new  problem  is  solved  by  the  conventional  exterior 
penalty  function  method,  because  this  is  believed  to  be  one 
of  the  most  efficient  algorithms  for  the  solution  of  such 
equality  constrained  problems  [35] . 


II.  HEAT  EXCHANGER  ANALYSIS 


A.  INTRODUCTION 

To  meet  the  objectives  of  this  thesis,  an  analysis  program 
for  an  air-cooled  heat  exchanger  must  be  coupled  with  a 
numerical  optimization  scheme  to  produce  a  complete,  detailed 
design  package.  COPES/CONMIN  has  greatly  simplified  this  task 

This  analysis  program  must  be  written  in  subroutine  form, 
titled  ANALIZ,  and  organized  into  three  segments:  input, 
execution,  and  output.  The  analysis  subroutine  must  also: 

1.  Take  into  account  the  variation  of  the  heat-transfer 
coefficients  and  differential  pressure  drop  with 
temperature  and/or  length  of  flow  path. 

2.  Be  iterative  free,  if  possible. 

3.  Be  written  in  such  a  manner  that  the  optimizer  will 
play  a  role  in  surface  selection. 

With  the  number  of  design  variables  approaching  the  prac¬ 
tical  limit,  the  importance  of  an  iterative  free  analysis 
subroutine  cannot  be  over-emphasized.  The  reason  being,  that 
at  the  beginning  of  each  design  iteration  in  CONMIN,  the  calcu¬ 
lation  of  all  gradients  (each  design  variable  and  active  con¬ 
straint)  requires  a  complete  pass  through  ANALIZ.  Therefore, 
the  computational  time  required  by  ANALIZ  directly  affects  the 
time  required  to  reach  the  optimum. 

B .  PROBLEM  FORMULATION 

The  air-cooled  heat  exchange:  is  shown  in  Figure  10.  A 
cross-flow  arrangement  with  both  fluids  unmixed  was  chosen. 
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Cool  air  enters  the  heat  exchanger  at  temperature  Tc 

pressure  p  ,  and  constant  specific  heat  c  .  The  cool  air 

Pa 

makes  one  pass  through  the  exchanger  as  it  flows  over  an  isos¬ 
celes  pitched  bank  of  finned  tubes.  The  air  is  heated  by 

water,  in  single  phase,  at  an  entering  temperature  of  T,  , 

ni 

and  constant  specific  heat  c 

Pw 

The  analysis  of  the  air-cooled ,  cross-flow  heat  exchanger 
centers  about  the  first  law  of  thermodynamics  and  on  the  heajt 
transfer  equation.  These  equations  as  they  apply  to  the 
exchanger  of  Figure  10  are  as  follows: 


macn  <TC  -Tc  ) 
a  pa  c2  C1 


(12) 


me  <T.  -T.  > 

W  Pw  hl  h2 


(13) 


Qc  =  U  A  AT 
5  mm 


(14) 


where  m  ,  T  and  m  ,  T.  are  the  fluid  mass  flow  rates  and 

a  c  2  w  “2 

exit  temperatures  of  air  and  water  respectively.  Um  is  the 
true  mean  overall  heat  transfer  coefficient  based  on  the  out¬ 
side  root  tube  area,  A  is  the  total  heat  transfer  surface 

area  of  the  exchanger  used  to  compute  U  ,  and  AT  is  the 

mm 

mean  temperature  difference  of  the  given  exchanger. 

The  object  of  the  analysis,  therefore,  is  to  determine 
•  •  • 

Q3»  Q,j  and  Q5  given  an  initial  listing  of  values  for  the 
design  parameters.  The  list  includes  the  following: 


Tubeside  mass  flow  rate,  lbm/hr 

Entrance  temperature  of  hot  stream,  °F 

Exit  temperature  of  hot  stream,  °F 

Specific  heat  of  hot  fluid,  BTU/lbm  -  °F 

Air  mass  flow  rate,  lbm/hr 

Entrance  air  temperature,  °F 

Exit  air  temperature,  °F 

Entrance  air  pressure,  psi 

Specific  heat  of  air,  BTU/lbm-°F 

Tube  inside  diameter  (ID) ,  in. 

Tube  outside  diameter  (OD) ,  in. 

Fin  height,  in. 

Fin  thickness,  in. 

Fin  spacing,  in. 

Transverse  Pitch,  in. 

Longitudinal  Pitch,  in. 

Bank  height,  in. 

Bank  width,  in. 

Cross-flow  arrangement 
Fin  type 
Number  of  rows 
Number  of  passes 

Given  heat  transfer  rate,  BTU/hr 

Thermal  Conductivity  of  tube  material,  BTU/hr .-ft.- °F 
Thermal  Conductivity  of  fin  material,  BTU/hr. -f t .-°F 
Among  the  initial  listing  of  design  parameters  above, 
there  are  parameters  that  are  known  and  will  remain  constant 


throughout  the  design  problem.  Also,  there  will  be  those 
parameters  that  are  unknown  and  can  vary,  i.e.  design  varia¬ 
bles.  On  the  way  to  determining  the  various  heat  transfer 
rates,  other  information  will  have  been  computed.  This  inform¬ 
ation  includes  the  objective  function,  constraining  functions, 
and  other  design  data,  such  as  the  number  of  tubes  per  verti¬ 
cal  row. 

The  optimizer  will  then  manipulate  the  design  variables 
in  order  to  find  an  optimum,  while  at  the  same  time,  performing 
a  heat  balance,  that  is: 

Q  =  Q3  =  Q4  =  Q5  (15) 

where  Q  may  be  some  given  heat  transfer  rate. 

C.  PERFORMANCE  CALCULATION  PROCEDURE 

With  the  temperatures,  mass  flow  rates,  and  specific 
heats  all  specified  in  the  listing  of  design  parameters; 
whether  they  be  constant  or  variable,  the  only  unknown  quan¬ 
tities  on  the  right  hand  side  of  equations (12)  through  (14), 

are  U  ,  A,  and  AT  .  They  will  be  determined  as  shown  in 
m'  m  J 

Figure  1. 

1.  Mean  Temperature  Difference  (MTD) 

For  many  flow  arrangements,  various  approahces  for 
determining  MTD,  mainly  using  diagrams,  are  available  [36], 
which  have  proven  very  useful  in  manual  design  efforts.  For 
computerized  design,  however,  an  explicit,  approximate  equation 
is  desirable  in  order  to  achieve  a  fast,  sufficiently  accurate 
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calculation  of  the  mean  temperature  difference  of  a  given 
flow  arrangement. 

Roetzel,  et  al.  [37],  presented  such  an  approximate 
equation  together  with  empirical  coefficients  for  nine 
counter-current  cross-flow  arrangements  as  they  apply  to 
air-cooled  heat  exchangers. 

Roetzel  used  the  familiar  equation  for  the  MTD  of  the 
given  flow  arrangement  AT : 


AT  =  F  •  AT. 
rn  l 


where  AT,  is  the  limiting  case  of  pure  countercurrent  flow: 

<'m 

<T  -T  )  -  <Th  -T  > 

’  -3--i - . -*-■  (17 ) 

m  (T.  -T  ) 

1n  hl  2 

1  (T.  -  T  ) 

h2  °1 

and  F  is  a  correction  factor  determined  by  a  different  set 
of  coefficients  for  each  flow  arrangement.  Roetzel  reported 
the  following  function  suitable  for  F: 


*11  ** 

=  1  -  l  l  ®i  kU-Vi 

i=l  k=l  1,K  r 


sin(2iarctan  R)  (18) 


where  is  the  dimensionless  LMTD , 

m 


T  _  T 
hl  C1 
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and  the  coefficients  a.  .  of  the  approximating  equation  (18), 

1  ,  K 

were  calculated  using  a  standard  least  squares  estimation 
program  [38]  ,  and  are  reported  in  reference  [37]  .  The  assump¬ 
tion  that  both  streams  were  unmixed  was  used  in  their  calcu¬ 
lation  . 

When  more  than  four  tubeside  passes  are  used,  it  is 
assumed  that  the  heat  exchanger  has  approached  the  limiting 
case  of  pure  counterflow  and  F  is  set  equal  to  one  [2] . 

Having  determined  the  MTD,  the  remainder  of  the  analy¬ 
sis  procedure  follows  Roetzel's  [39]  general  approximation 
method  for  determining  the  mean  overall  heat  transfer  coeffi¬ 
cient,  Um,  for  any  flow  arrangement  while  taking  into  account 
the  variation  of  the  heat  transfer  coefficients  and  the 
pressure  drop  with  temperature  and/or  length  of  path. 

Before  continuing  with  specific  analysis  procedures, 
a  brief  summary  of  Roetzel's  general  approximation  method  is 
in  order . 

The  local  overall  heat  transfer  coefficient  based  on 
the  outside  root  tube  area  can  be  written  as  follows: 


U  = 


A~  I  A  in (r/r . )  \ 

ol,o  o  l  ,  1 

A.  H.  2 irkL  H  n* 

11  of 


(21) 


where  is  the  efficiency  of  the  extended  surface. 

In  order  to  determine  the  individual  convection  heat 


transfer  (film)  coefficients,  and  Hq,  according  to  the 
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conventional  methods,  the  coefficients  would  be  considered 
constant,  and  the  necessary  fluid  properties  for  their  calcu¬ 
lation  would  be  evaluated  at  some  mean  bulk  temperatures, 

T,  and  T 
hb  cb 

However,  the  film  coefficients  are  not  constant,  but 

vary  with  temperature  and/or  length  of  flow  path.  Roetzel 

has  taken  these  variations  into  account  with  the  use  of 

corrected  reference  temperatures.  Two  sets  of  corrected 

reference  temperatures  are  determined:  T,  ,  T  and 

nI  CI 

T.  ,  T  .  Therefore,  for  each  set  of  corrected  reference 
nII  CII 

temperatures,  the  film  coefficients  are  determined  in  the 
conventional  manner  using  the  reference  temperatures  in  place 
of  the  bulk  temperatures . 

With  the  film  coefficients,  H.  ,  H  ,  H.  and  H  , 

X1  °I  XII  °I  I 

two  local  overall  heat  transfer  coefficients  can  be  calcu¬ 
lated  from  equation  (21),  and  U 

Finally,  the  true  mean  overall  heat  transfer  coeffi¬ 
cient  is  calculated  as: 


JL  -  i[_L  +  __L_ 

um  2  01  Un 


(22) 


2 .  Determination  of  Reference  Temperatures 

The  reference  temperatures  for  a  pure  counterflow 
heat  exchanger  must  first  be  determined  from  [39] : 


Th. 


=  T.  +  (T.  -T.  ) 

“2  hl  h2 


&Tt  -  (T,  -T  ) 

1  h2  C1 
(T.  -T  )  -  (T.  -T  ) 
hl  C  2  h2  C1 


(23) 
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4TII  -  <Th '%  > 

ThII  "  Th,  +  <Th1'Th,)  (T.  -T  )  -  (Th  -T  ; 

II  2  12  ”l  c2  “2  C1 


ATI  "  (Th  -Tc  ) 

Tc  =  Tc.  +  (Tc  ’V  (Th  ~Tc  7  "  (^h  ~Tc  } 
X  X  2  1  c  2  rij 


T  +  (T  -  T  ) 
C1  c2  C1 


AT__  -  (T.  -  T  ) 

II  xl  2  C  ^ 

(T.  -T  )  -  (T.  ~T  ) 
hl  C2  h2  C1 


where 


AT  =  (T  -T  )*78868  .  (T  -t  )*21132  (26a) 

1  hl  c2  h2  C1 


AT  =  (T  -T  )‘21132  •  (T  -T  )‘78868  (26b) 

11  hl  c2  h2  C1 


Equations  (26a)  and  (26b)  were  derived  by  Roetzel  in  references 
[40]  and  [41]  .  For  the  special  case  where  the  fluid  heat 
capacity  rates,  C,  are  equal,  that  is. 


m,  C  =  C 


h  p 


C  P  r 


the  term 


AT.  -  (T,  -T  ) 

3  h2  C1 

[Th  "Tc  >  '  (Th  -Tc  } 
n2  c2  n2  c1 


of  equations  (23)  through  (26),  with  j  =  I  and  II,  becomes 


.78868  or  .21132  respectively. 
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Correction  of  Reference  Temperatures  for  Given 


With  inlet  and  outlet  temperatures  fixed,  pure  counter¬ 
flow  yields  the  highest  mean  temperature  difference.  There¬ 
fore,  for  any  other  arrangement,  the  temperature  difference 
would  be  smaller.  Thus,  the  corrections  are  applied  in  the 
following  manner  [39] : 


T  =  T'  +  4/ 

CII  CII  CII 

where  +,  and  4/  are  the  temperature  corrections  of  the 
*  ■:  k 

hot  and  cold  streams,  respectively.  The  corrections  are 
calculated  as  follows: 


The  corrected  reference  temperatures  are  now  used  to 
determine  the  thermal  conductivities  and  absolute  viscosities  of 
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fluids  for  later  use  in  the  calculation  of  the  film  coeffi¬ 
cients.  Thermal  conductivity  and  viscosity  data  are  usually 
presented  in  tabular  form.  However,  for  use  on  the  computer, 
an  explicit,  simple  approximate  equation  with  temperature  as 
the  independent  variable,  was  desirable.  Water  and  air  were 
chosen  as  two  fluids  that  were  likely  to  be  involved  in  an 
air-cooled  heat  exchanger  design.  Figures  11,  12  and  13 
indicate  that  the  thermal  conductivities  of  air  and  water 
and  the  viscosity  of  air  can  be  approximated  by  a  second  order 
polynomial.  The  viscosity  of  the  hot  tubeside  fluid,  water, 
must  be  treated  specially,  due  to  the  following  considerations. 

Calculations  for  the  tubeside  heat  transfer  coeffi¬ 


cient  in  the  laminar,  transition  and  turbulent  regions  all 

14 

involve  the  correction,  {y/ywall)*  [ 42].  In  the  past,  an 
iterative  procedure  was  required  to  determine  the  inside  tube 


wall  temperature  in  order  to  evaluate  uwaly  Roetzel  [43], 
has  developed  an  iteration-free  method  for  determining  this 


correction.  Roetzel's  method  assumes  that  the  tubeside 


fluid's  viscosity  follows  Andrade's  Correlation,  that  is: 


v 


a 


e0/T 


Using  the  viscosity  data  of  saturated  water  from  reference 
[44] ,  the  coefficients  were  determined  through  regression 
analysis  with  the  resulting  equation: 
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where  temperature  is  in  degrees  Rankine  and  viscosity  is  in 
lbm/ft-hr.  Equation  (27)  was  plotted  next  to  the  experimental 
data  with  results  that  indicate  that  water  does  indeed  satisfy 
Andrade's  equation,  see  Figure  14. 

4 .  Uncorrected  Tubeside  Heat  Transfer  Coefficient 
The  tubeside  heat  transfer  coefficient  is  calcu¬ 
lated  from  one  of  three  Nusselt-type  empirical  equations 
[42],  as  follows : 

For  laminar  flows,  Reynolds  Number  <  2,100 


1.86  Re1''3  Pr1/3  &)  ^  <-l_,  -14 

L  ywall 


(28) 


For  transition  regions,  2,100  <  Re  <  10,000 


[Re 


2/3 


-  125] [1  + 


£>2/3 


]  Pr 


1/3 


-) 


.14 


wall 


(29) 


For  turbulent  flow.  Re  >  10,000 


HD. 

i 

k 


.023  Re'8  Pr 


.1/3 


(7") 

w 


,14 


(30) 


Therefore,  before  any  calculations  can  even  begin, 
the  Reynolds  Number,  Re,  must  be  computed  to  determine  the 
type  of  flow. 

The  Reynolds  Number  will  be  calculated  as: 


D . 

l 
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is  supplied  in  the  initial  list  of  parameters  and 
u ^  is  calculated  from  equation  (27)  .  This  leaves  only  the 
cross-sectional  flow  area,  Ax»  to  be  determined,  where: 


The  number  of  tubes  per  pass,  N^,  is  determined  geometrically, 
having  been  given  the  tube  bank  height,  number  of  rows,  num¬ 
ber  of  passes,  transverse  pitch,  diameters,  and  fin  height 
initially.  Figure  15  details  the  procedure  of  finding  N  . 
Mathematically , 


N 


P 


[h  .  Df 

pass  1 


and  would  be  a  rational  number.  This  rational  number  is 
then  truncated  to  an  integer  number  of  tubes . 

With  Re  calculated,  the  Prandtl  Number, 


is  computed  from  given  and  previously  determined  thermo-physical 
properties . 

The  uncorrected  tubeside  heat  transfer  coefficients, 

14 

H!  ,  that  is;  H.  ,  without  the  factor  (y/u  ) *  ,  can  be 

calculated  from  the  proper  choice  of  equations  (28)  through 
(29)  . 
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5.  Calculation  of  Wall  and  Associated  Resistances 


as : 


Equation  (21)  can  be  written  in  a  more  general  form 


U 


R.  +  R  , ,  "+  R 
i  wall  o 


(31) 


where  R^,  Rw  and  Rq  are  the  inside  wall  and  outside  heat 
transfer  resistances  respectively.  Additional  resistances, 
such  as  contact  or  fouling,  can  also  be  added  here. 

In  comparing  equations  (21)  and  (31)  ,  the  resistances 
can  be  computed  as : 


R. 


l 


j 


Rwall 


D.  H. 
i  i 


j 

ln(Do/Di) 


^wall 


Ola) 


(31b) 


j 


(31c) 


6 .  Calculation  of  Airside  Heat  Transfer  Coefficients 
In  order  to  have  the  optimization  program  play  a 
significant  role  in  the  selection  of  an  optimized  surface 
for  a  finned  tube  heat  exchanger  (which  is  one  of  the  objec¬ 
tives  of  this  project) ,  an  explicit  equation  for  Hq  involving 
tube  and  bank  design  parameters  as  independent  variables  is 
a  necessity. 

In  the  past,  comparison  methods,  as  described  by  Shah 
[45] ,  were  used  to  choose  the  best  surface  from  a  list  of 
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surfaces  for  which  experimental  heat  transfer  and  friction 

data  existed.  The  data  is  presented  in  graphical  form, 

where  Friction  and  Colburn  Factors  are  plotted  as  a  function 

• 

of  Reynolds  Number.  With  J,  m  ,  c  and  Pr  known,  the  film 

a  p 

coefficient  can  be  computed. 

Therefore,  in  previous  heat  exchanger  optimization 
programs,  a  given  surface  described  by  its  pitch,  outside 
fin  diameter,  fin  thickness,  fin  spacing,  and  outside  tube 
diameter,  had  to  be  chosen  bef orehand .  After  the  surface 
configuration  had  been  selected,  expressions  for  f  and  J 
were  obtained  by  fitting  polynomials  to  the  experimental 
data  described  earlier. 

Briggs  and  Young  [47]  provide  a  means  of  getting 
past  this  obstacle  with  an  improved  convection  heat  transfer 
correlation  for  air  flowing  across  triangular  pitch  banks  of 
high  finned  tubes. 

Briggs  and  Young  expanded  the  work  of  Ward  and  Young 
[48],  to  cover  a  total  of  18  differently  configured  finned 
tube  banks  in  order  to  determine  the  effect  of  fin  thickness 
and  tube  pitch  on  the  airside  heat  transfer  coefficient.  The 
heat  transfer  data  for  the  high-fin  tube  banks  were  correlated 
to  gives 

Nu  =  — -  =  .1378  Re*718  Pr1/3  (f-)'296  (32) 

j 


where  s  is  the  distance  between  adjacent  fins  and  i  is  the  fin 
height.  Equation  (32)  is  based  on  tubes  having  a  wide  range 


of  fin  heights,  fin  thicknesses,  fin  spacing  and  outside 

tube  diameter  and  can  be  used  to  predict  H  for  a  bank  of 

°j 

tubes  six  rows  deep.  Figure  11  of  reference  [48]  is  used  to 
correct  H'  for  banks  of  other  than  six  rows. 

°j 

7 .  Calculation  of  Fin  and  Surface  Efficiencies 

In  order  to  calculate  the  outside  heat  transfer 
resistance,  which  will  be  used  to  calculate  the  correction 
for  the  tubeside  heat  transfer  coefficient  and  finally,  the 
local  overall  heat  transfer  coefficients,  the  extended  sur¬ 
face  efficiency,  nf,  must  be  computed. 

The  surface  efficiency  accounts  for  the  temperature 
drop  from  the  root  to  the  tip  of  the  fin,  due  to  the  thermal 
resistance  of  the  fin  material.  Thus,  even  though  the  heat 
transfer  has  been  increased  by  the  additional  area  of  the 
extended  surface,  the  area  is  not  as  effective  as  if  it  were 
to  be  all  at  the  root  temperature . 

The  surface  efficiency  can  be  expressed  as  [49] : 

^f 

nf  =  1  •  X  (1_1>) 

where  <? ,  the  fin  efficiency  for  a  radial  fin,  can  be  found 
from  [49] : 


2 


r 

o 


I.  (mr  _)  k.  (mr  )  -  k.  (mr, )  I.  (mr  ) 

r  1  f  1  O _ 1  f  1  O  -I 

lIQ (mro)k1 (mrf )  +  1^ (mrf ) kQ (mrQ) 1 


The  finned  area,  A^,  and  the  total  heat  transfer  area.  A,  are 
computed  as  follows: 
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As  »  Nt  w  nf  i(Df2  -  Do2) 


A  =  NTw[^(Df2 -Dq2)  +  (1  -  nft)DQiT3 


With  Rq  calculated  from  equation  (31c) ,  the  correction  to 
the  tubeside  film  coefficient  can  be  made. 

8 .  Correction  of  the  Tubeside  Heat  Transfer  Coefficient 
With  the  tubeside  heat  transfer  coefficient  dependent 
on  the  wall  temperature,  the  dependence  has  either  been 
neglected,  or  the  wall  temperature  has  been  calculated  with 
an  iterative  technique  in  the  past.  Roetzel  [43]  ,  has  pro¬ 
posed  an  improved  iterative-free  method  for  finding  the  "Seider- 

14 

Tate"  correction,  (u/uwall) '  ,  when  the  tubeside  fluid 

obeys  Andrade's  viscosity  correlation. 

From  Roetzel's  formulation: 


14 

wwall  j 


u .  u  . 

-  -i  *  <-4~  +  vj> 


1/2 


where 


u . 
3 


.007  In ( 


—  j)^- 


3 


B 


1  -  Tc/Th. 
_ —J - ^ 

Th./Twall *  1 


Vh, 


-  1 


v . 
3 


wwall .  x 

■°7 

3  3 


1  -  T  /T 

hi 

Tv  /T* -1 


h 

3 


wall 


B:Th. 


UL 


45 


where 


B  . 
3 


A  (R  . , 
wall 


+  V 

J 


(Note  that  all  temperatures  are  in  degrees  absolute.) 

All  the  parameters  on  the  right-hand  side  of  equation 
(21)  are  now  available.  Therefore,  the  two  reference  over¬ 
all  heat  transfer  coefficients,  Uj  and  u^,  can  be  calculated 
The  mean  overall  heat  transfer  coefficient  follows  easily 
from  equation  (22) . 

With  U^,  Q5  can  be  calculated,  with  the  heat  balance 
to  be  performed  by  the  optimizer. 

9 .  Calculation  of  Air  and  Tubeside  Pressure  Drops 

The  final  calculations  before  computing  objective  and 
constraint  functions  for  the  optimization  problem  involve 
the  pressure  drops  in  the  heat  exchanger. 

The  basic  equations  that  will  be  used  for  the  calcu¬ 
lation  of  the  pressure  losses  are  as  follows: 
a.  Tubeside  [42] 


Ap 


•  2 

f.  m,  w  n 
i  h 


1  2  9C  Ax  p  Di  (IT) 

c  XT  1  ww 


z  + 


(n  -  1)  m  2 
-_P  .  ^  h 


(33) 


2  gc  V  p 


where  z  *  .14  below  Re  =  2100  and  z  =  .25  for  Re  greater  than 

2100. 

b.  Air side  [50] 


*Pr 


•  2 

f  n  m  z 
ore 
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(34) 
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where  A^f  is  the  minimum  flow  area. 


The  friction  factor  for  the  tubeside  flow,  f^,  is 
taken  from  Figure  9.5  of  reference  [42],  assuming  fully 
developed  flow.  For  use  on  the  computer,  an  explicit  expression 
for  f .  was  obtained  by  fitting  a  line  and  an  exponential  to 
the  experimental  data  of  Figure  9.5.  This  follows: 


Re  <_  1000  fi  »  .5/Re 

Re  >  1000  f.  =  .003215  Re"’2694 

Just  as  in  the  case  of  the  airside  film  coefficient, 

for  surface  optimization  on  the  computer,  an  explicit  equation 

for  the  airside  friction  factor,  f  ,  is  a  necessity.  Robinson 

o 

and  Briggs  [50]  ,  presented  such  an  expression  for  f  for  air 
flowing  across  triangular  pitch  banks  of  finned  tubes.  Robinson 
and  Briggs'  work  closely  parallels  that  of  Briggs  and  Young 
[47] .  The  Robinson-Briggs  Correlation: 


f 


o 


18.93  Re"*316 


(34a) 


covers  the  range  of  tube  sizes  and  pitches  used  in  air-cooled 
heat  exchangers  [50] . 

Therefore,  with  Rej  and  Re^,  the  four  peference 
pressure  drops,  Ap.  ,  and  Ap  ,  may  be  computed  from  equations 

1j  °j 

(33)  and  (34).  Following  Roetzel's  general  approximation 


method: 
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For  a  gas,  an  additional  correction  is  needed  because  the 
density  in  equation  (33)  is  strongly  dependent  on  pressure, 
which  is  changing  through  the  exchanger.  Using  the  inlet 
pressure  as  reference: 
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All  the  necessary  information  from  an  analysis  view¬ 
point  has  now  been  calculated.  Functions  needed  for  the 
numerical  optimization  process  shall  follow. 

10 .  Objective  and  Constraint  Functions 

The  objective  functions  available  for  minimization 
are  defined  as  follows: 

a.  Volume  =  wh[D,  +  (n  -1)P_  cos  3J 

r  r  L 

where 


©  =  arcsin  (P  /2P_ ) 

t  Li 


(36) 


b. 


(Df2  "  Do2)  n  t)D  -1 

Heat  Transfer  Area  =  NTw[ - g -  +  o  1 


c.  Air  Horsepower  = 


^P. 


m 


48 


d.  Airside  Pressure  Drop 

e.  Tubeside  Pressure  Drop 

Limitations  were  imposed  on  the  following  inequality 
constraints  in  order  to  keep  the  design  within  practical 
physical  bounds: 

a.  The  diameter  ratio, 

Df 

D RATIO  =  =p 
O 

must  be  kept  reasonable.  This  can  also  be  accomplished 
to  some  extent  by  placing  side  constraints  on  the 
design  variables,  i  and  Dq.  See  Figure  16. 

b.  The  optimizer  must  be  prevented  from  driving  the 
tube  thickness. 


TUBTH  *  (DQ-Di)/2 


to  zero. 

c.  The  tubes  must  be  kept  from  touching  in  both  the 
longitudinal  and  transverse  directions. 


TOUCHN  =  Df  -  Pfc 
TOUCHL  =  Df  -  PL 

TOUCHL  and  TOUCHN,  therefore,  must  be  kept  below  zero. 


d.  Reasonable  temperature  profiles  must  be  maintained 
at  both  ends  of  the  exchanger.  See  Figure  10-7  of 
Ref,  [44], 


PROFH 


1 


PROFC 


that  is,  PROFH  and  PROFC,  must  be  negative. 
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e.  The  number  of  tubes  per  vertical  row,  VROWR,  shall 
not  be  allowed  to  go  below  2 . 

f.  The  free  face  area,  that  is,  the  minimum  flow  area 
for  air,  must  obviously  be  greater  than  zero. 


DELSFF  =  [projected  tube  area]  -  hw 


where  the  projected  tube  area,  STOTAL,  is 

D,  tw 

STOTAL  «  N  [D  w  +  - - ] 

r  o  s 

and  N  is  the  number  of  tubes  per  vertical  row.  DELSFF 
must  fie  less  than  zero. 

g.  The  airside  and  tubeside  pressure  drops  must  be  kept 
within  design  constraints, 

h.  From  experience,  the  argument  of  the  square  root  of 
equation  (35) , 

ARG5  =  1  - 

P1 

has  been  driven  below  zero.  It  must  therefore  be 
constrained . 

i.  To  maintain  an  isoscele  pitch  bank,  the  angle,  3, 

as  shown  in  Figure  17  and  as  defined  in  eq.  (36)  ,  must 
be  constrained.  One  such  constraining  value  is: 


0  =  arcos  (D,/2P  )  (37) 

m  r  L 

The  arguments  of  the  arcs in  and  arcos  of  eq.  (36)  and 
(37)  must  be  constrained  from  going  beyond  1. 

j.  The  heat  balance  described  by  eq.  (15)  is  performed 
by  three  equality  constraints: 


Q RATIO  =  Q5/Q 

(38) 

QRATIOl  =  Q3/Q5 

(39) 
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QRATI02 


(40) 


-  VQ3 

All  constraints  are  set  equal  to  1. 

Besides  the  constraints  described  in  paragraphs  a.  through 
j.,  above,  lower  and  upper  bounds  are  placed  on  the  design 
variables  to  assure  a  reasonable  design. 


IV.  RESULTS 


A.  BACKGROUND 

Case  studies  were  chosen  as  the  best  way  to  test  the 
capabilities  of  the  program  for  Heat  Exchanger  Design  using 
Numerical  Ogtimization  (HEDSUP) .  The  design  problems  posed 
were  made  as  realistic  as  possible. 

1 .  Capabilities 

HEDSUP  currently  has  the  capability  to  design  for 
nine  different  configurations  of  triangular  pitch  banks  of 
finned  tubes: 


TYPE 

1 

— 

1 

ROW, 

1 

PASS 

TYPE 

2 

- 

2 

ROW, 

1 

PASS 

TYPE 

3 

- 

3 

ROW, 

1 

PASS 

TYPE 

4 

- 

4 

ROW, 

1 

PASS 

TYPE 

5 

- 

2 

ROW, 

2 

PASS 

TYPE 

6 

- 

3 

ROW, 

3 

PASS 

TYPE 

7 

- 

4 

ROW, 

2 

PASS 

TYPE 

8 

- 

4 

ROW, 

2 

PASS 

TYPE 

10 

_ 

PURE  COUNTERFLOW 

TYPE  10  will  include  exchangers  with  a  configuration  of  n  rows, 
n  passes,  where  n  can  go  from  five  to  20. 

The  banks  must  be  constructed  of  high-finned  tubes 
( l  >  .0625  in.  [42])  with  the  fins  having  a  rectangular  pro¬ 
file  of  constant  thickness.  Additional  profiles  can  be  inserted 
into  HEDSUP  quite  simply,  provided  that  its  fin  efficiency  can 
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be  expressed  explicitly  as: 


n  =  f(2,t,H,k) 

See  Subroutine  FINEFF  of  the  program  listing.  Appendix  D. 

At  present,  HEDSUP  can  provide  the  design  parameters 
for  an  air-cooled  heat  exchanger  optimized  for  any  one  of  the 
following  design  objectives: 

(a.)  Minimum  Volume 

(b.)  Minimum  Heat  Transfer  Surface  Area 
(c.)  Minimum  Air  Horsepower 
(d.)  Minimum  Air side  Pressure  Drop 
(e.)  Minimum  Tubeside  Pressure  Drop 
Additional  design  objectives  can  be  used,  provided  that  they 
can  be  expressed  explicitly  as  a  function  of  the  design  varia¬ 
bles  and  they  are  added  to  the  common  block.  It  should  also 
be  pointed  out  that  any  design  variable  may  simultaneously  be 
a  design  objective  as  long  as  it  conforms  to  the  restrictions 
of  both.  For  example,  an  exchanger  may  be  designed  for  mini¬ 
mum  bank  height. 

The  airside  fluid  is  restricted  to  dry  air.  The  tube- 
side  fluid  is  presently  limited  to  water  in  single  phase. 

Other  tubeside  fluids  can  be  used  by  HEDSUP,  provided  that 
their  visocisites  obey  Andrade's  Law  and  the  fluid's  thermal 
conductivities  can  be  expressed  explicitly  as  a  function  of 
temperature.  The  fluid's  specific  gravity  would  also  have  to 
be  placed  in  the  denominator  of  eq.  (33)  . 
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B.  CASE  STUDIES 


1 .  Case  One 

a.  Problem  Formulation 

An  air-cooled  heat  exchanger  is  to  be  designed 
for  minimum  volume  with  a  heat  transfer  rate  of  10,000,000 
BTU's  per  hour.  Water  is  to  be  coaled  from  200 °F  to  125 °F. 
Dry  air  will  enter  the  exchanger  at  95°F  and  leave  at  130°F 
Specifications  call  for  a  fan  that  can  produce  a  pressure 
difference  of  two  inches  of  water. 

b.  Design  Variable  Framework 

From  the  list  of  design  parameters  in  Section 
III.B,  the  design  variables  must  be  singled  out,  identified 
to  COPES,  and  given  side  constraints.  All  parameters  must 
be  given  an  initial  value.  Only  the  values  of  the  design 
variables  will  change. 

Assuming  constant  specific  heats. 


c  =1.0  BTU/lbm-°F 
Pw 


c  =  .24  BTU/lbm-°F 

P* 


the  mass  flow  rates  of  both  fluids  can  be  determined  from 
eqs.  (12)  and  (13)  ,  because  the  heat  transfer  rate  and  tem¬ 
perature  differences  are  given, 


m. 


w 


m 


a 


133333.  lbm/hr 


190476.2  lbm/hr 


A  cross-flow  arrangement,  fin  profile,  tube 
material,  and  fin  material  must  be  chosen. 


The  design  variables  for  this  example  are,  there¬ 


fore? 


,232  <  Di  <  2.325  in. 


,24  <  D  <2.5  in. 

o 


,0625  in.  <  i  <  <*> 


,01  <  t  <  .0235  in. 


,08  <  S  <  .125  in. 


0.0  <  P_  <  4.0  in. 

h 


0.0  <  <  4.0  in. 


=  2.0  inches 


=  2.5  inches 


,46  inches 


=  .023  inches 


=  .111  inches 


=  2.125  inches 


4 . 0  inches 


0.0  <  w  <  500  in. 


=  490  inches 


0.0  <  h  <  500  in. 


=  350  inches 


The  side  constraints  on  the  design  variables  are  of  a  practi¬ 
cal  nature  with  the  exception  of  the  lower  bounds  on  fin 
height.  Recall  that  the  use  of  eq.  (32)  is  restricted  to 
high-fins.  High-fins  will  also  tend  to  keep  the  fluid  unmixed, 
which  was  an  assumption  used  when  defining  the  coefficients 
a.  t,  used  in  eq.  (18) . 

c.  Constraint  Framework 

From  the  problem  statement,  the  airside  pressure 
drop  must  be  less  than  two  inches  of  water  or  .0722  psi. 
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0  <  Apa  <  .0722  psi 
From  a  practical  standpoint: 
0.0  <  Q  <  1.3 

1.0  <  DRATIO  <  2.5 


.018  <  TUBTH  <  .18  inches 


where 


-oo  <  p  <  .14  psi 

rw  c 

The  equality  constraint 


Q 


Qq  =  U  A  AT 
j  m 

Q  =  10,000,000  Btu/hr 


will  satisfy  the  heat  balance. 

The  bounds  on  TOUCHN,  TOUCHL  ARG5,  ARG7  ,  ARG8 , 
DELSFF  and  VROWR  were  discussed  in  Section  I II. C. 10.  Con¬ 
straints  on  PROFH  and  PROFC  were  unnecessary  because  all 
temperatures  were  specified  in  the  problem  statement. 

As  with  design  variables,  constraints  must  be 
identified  to  COPES  by  location  in  the  common  block.  See 
the  global  catalog.  Appendix  A. 
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d.  Methodology 


Ideally,  a  three  dimensional  design  matrix  can 
now  be  constructed  of  optimum  exchanger  designs  with  minimum 
volumes.  The  matrix  would  be  constructed  by  first  holding 
the  tube  and  fin  materials  constant  and  varying  the  configura¬ 
tion,  i.e.  Type  1,  Type  2,  Type  3,  etc.  Next  the  tube  material 
would  be  varied  with  the  fin  material  and  exchanger  configura¬ 
tion  held  constant  and  so  forth.  However,  for  case  study  one, 
the  tube  material  will  be  chosen  as  copper,  k  =  200  BTU/f t-hr-°F , 
and  the  fin  material  will  be  aluminum,  k  =  118  BTU/ft-hr-°F. 

Also,  in  order  to  simulate  an  actual  trade  off 
study,  the  constraint  framework  will  be  fixed  throughout  the 
individual  case  studies. 

Problems  arise  in  constructing  the  matrix  when 
trying  to  determine  the  true  minimum  volume  design  for  each 
configuration.  Unfortunately,  the  choice  of  initial  design 
parameters  (starting  point) ,  coupled  with  the  input  parameters 
for  ALMM,  will  sometimes  lead  to  entirely  different  optimum 
designs  with  volumes  differing  by  over  100%. 

The  ALMM  parameters  include  the  initial  multiplier, 
CC,  the  multiplication  factor,  CMULT,  and  the  maximum  multi¬ 
plier  value,  CCMAX.  Experience  has  shown  that  setting 

CMULT  =  2.0 

CCMAX  =  1000. 

will  suffice  for  almost  all  problems.  However,  there  is  much 
"artwork"  involved  with  the  choice  of  CC.  From  experience. 


* 
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an  initial  multiplier  of  10  works  well  when  starting  far 
from  the  equality  constraint,  i.e.: 

.8  <  Q RATIO  <  1.2 

However,  when  approaching  the  heat  balance,  i.e.,  the  equality 
constraint,  a  CC  =  100  is  necessary  for  the  heat  balance  to 
converge.  Equation  (15)  is  considered  satisfied  when  the  heat 
transfer  ratios  (eqs.  (38)  through  (40))  are  less  than  0.1 
percent. 

Therefore,  it  is  obvious  why  the  chosen  initial 
design  is  so  critical.  Together  with  the  choice  of  CC  it 
will  determine  how  and  to  what  design  the  optimizer  will  con¬ 
verge.  As  an  example,  see  Table  1.  Notice  the  calculated 
heat  transfer  for  the  initial  design.  This  value  is  the  product 
of  U^,  A,  and  ATm  calculated  by  ANALIZ  using  the  initial  design 
parameters,  some  of  which  are  mere  estimates.  Recall  it  is 
the  job  of  the  optimizer  to  vary  the  design  variables  in  order 

•  at  • 

to  bring  Q5  equal  to  Q,  and  Q4  and  at  the  same  time  minimize 
the  objective  function. 

In  order  to  remove  some  of  the  "artwork"  and 
try  to  ensure  a  true  optimum  design,  the  following  procedure 
is  recommended  to  generate  the  design  matrix: 

(1)  Begin  with  a  Type  2  configuration;  input  the  initial 
design  and  constraint  values  enumerated  in  sections 
(a)  through  (d)  above;  let  CC  =  10.;  execute. 

(2)  If  the  heat  balance  of  the  resulting  design  has  not 
converged,  but  is  within  20%,  use  the  design  results 
as  a  new  starting  point  and  set  CC  =  100  (if  the  heat 
balance  is  not  with  20%,  let  CC  =  10.). 
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TYPE  7  EXCHANGER  CONFIGURATION  OPTIMUM  DESIGN 


RUN 

I 

RUN  II 

DESIGN 

INITIAL 

OPTIMUM 

INITIAL 

OPTIMUM 

VARIABLES 

DESIGN 

DESIGN 

DESIGN 

DESIGN 

D^  ,  inches 

.  6786 

.  5392 

2.0 

1.168 

D  ,  inches 
o 

.7201 

.  5767 

2.5 

1.278 

l  ,  inches 

.1618 

.128 

.46 

.294 

t  ,  inches 

.0217 

.0208 

.023 

.023 

S  ,  inches 

.08 

.08 

.111 

.08 

Pfc  ,  inches 

1.044 

1.02 

4.0 

2.557 

PL  ,  inches 

1.044 

.841 

2.125 

1.899 

h  ,  inches 

342.1 

238.4 

350. 

172.9 

w  ,  inches 

286.6 

186.2 

490. 

253. 

Q  ,  BTU/hr 

21,988,624. 

9,999,925. 

20,650,352. 

9,999,811. 

Volume,  ft3 

213.12 

72.92 

553.2 

153.83 

Table  1 
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(3)  Repeat  step  (2)  until  convergence.  NOTE:  CC  may 
be  adjusted  up  to  150  when  approaching  convergence. 

The  design  should  converge  following  the  use  of 

CC  =  100.  If  too  much  adjustment  of  CC  is  necessary, 
reaching  the  optimum  from  the  starting  point  is 
unlikely. 

(4)  Ensure  a  minimum  design  by  beginning  step  (1)  with  a 
different  initial  design. 

(5)  After  finding  an  optimum  design  for  a  Type  2  con¬ 
figuration,  use  that  design  for  the  starting  point 
for  a  Type  3  configuration.  This  assures  a  reason¬ 
able  starting  point,  probably  close  to  the  optimum 
for  Type  3 . 

(6)  Repeat  steps  (1)  through  (3)  for  a  Type  3  configuration. 

(7)  Repeat  steps  (5)  and  (6)  for  the  remaining  configuration 

e.  Design  Matrix 

The  design  matrix  is  presented  in  Table  2.  The 
optimum  design  is  a  Type  Four  configuration  and  is  shown  in 
Figure  18.  Typically,  for  this  case  study,  when  starting  far 
from  the  final  design,  COPES/CONMIN  would  require  approxi¬ 
mately  1700  calls  to  ANALIZ  to  reach  an  optimum.  However 
when  beginning  from  a  reasonable  starting  point  with  CC  =  100, 
only  600  calls  were  needed.  Note  that  each  call  to  ANALIZ 


requires  approximately  .06  seconds  of  CPU  time  on  an  IBM  360/67 


2.  Case  Study  Two 

a.''-'-  Problem  Formulation 


An  air  heater  is  to  be  designed  to  fit  into  a 
space  8*  X  24'  X  4'.  The  heat  exchanger  is  to  heat  1,000,000 
lbm/hr  of  dry  air  from  75°F  to  130°F.  256,000  lbm/hr  of  water 

at  200 °F  is  available.  Design  the  heater  so  that  the  required 
air  horsepower  is  at  a  minimum. 
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b.  Design  Variable  Framework 

Assuming  constant  specific  heats. 


c  =1.0  BTU/lbm-°F 

?w 

c  =  .24  BTU/lbm-°F 

Pa 

eqs.  (12)  and  (13)  will  yield  the  required  heat  transfer 
rate  and  the  outlet  water  temperature. 

Q  *  m  C  (T  -  T  )  =  13.2  *  106  BTU/hr 

pa  c2  C1 

Cp  Thx  "  Q 

T.  =  - - - -  =  148.44  °F 


Therefore,  the  design  variables  are  as  follows: 


.232 

< 

D. 

1 

< 

2.325  in. 

as 

2.0  in. 

.25 

< 

D 

< 

2.5  in. 

D* 

2.5  in. 

o 

0 

.0625 

in.  < 

l 

<  00 

= 

.46  in. 

.01 

< 

t 

< 

.0235  in. 

t1 

= 

.023  in. 

00 

o 

• 

< 

S 

< 

.125  in. 

s1 

= 

.Ill  in. 

o 

• 

o 

< 

PL 

< 

m 

V 

= 

2.125  in 

o 

• 

o 

< 

Pt 

< 

CO 

pti 

* 

4.00  in. 
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c.  Constraint  Framework 


When  designing  for  minimum  horsepower,  the  opti¬ 
mizer  will  naturally  try  to  drive  the  design  to  a  maximum 
volume  in  order  to  reduce  air side  pressure  losses.  Therefore, 
it  is  reasonable  to  constrain  volume  as  follows: 

0  <  Volume  <  768  ft^ 

The  other  constraints  are : 

0.0  <  Q  <  1.3 

1.0  <  DRATIO  <  2.5 

.018  <  TUBTH  <  .18  in. 


0  <  Apw  <  .14  psi 

and  for  the  heat  balance: 


d.  Design  Matrix 

With  the  tube  and  fin  materials  fixed,  as  in  case 
one,  the  matrix  is  presented  in  Table  3.  The  optimum  design 
is  a  Type  4  configuration,  as  shown  in  Figure  19.  In  case 
study  two  when  starting  far  from  the  optimum,  the  optimizer 
called  ANALIZ  approximately  1900  times.  When  starting  close 
to  the  final  design,  for  example  using  the  design  for  a  Type  3 
configuration  as  a  starting  point  for  the  Type  4  design, 
COPES/CONMIN  only  required  596  calls. 


63 


Table  3.  CASE  II  DESIGN  MATRIX 


V.  CONCLUSIONS 


The  intent  of  this  investigation  was  to  couple  an  analysis 
program  with  a  numerical  optimization  scheme,  COPES /CONMIN, 
to  produce  a  complete,  detailed  design  program  for  an  air¬ 
cooled  heat  exchanger,  (HEDSUP) .  In  addition,  the  analysis 
program  was  to  be  written  such  that: 

1.  the  variation  in  the  film  coefficients  with  temperature/ 
length  of  flow  path  would  be  taken  into  account 

2.  the  surface  would  be  optimized 

3.  it  would  be  iterative  free  and  thus  minimize  the  CPU 
time  required  during  an  actual  trade  off  study. 

The  results  from  test  cases  using  ANALIZ  coupled  with 
COPES/CONMIN  in  its  present  form  were  unsatisfactory.  Although 
COPES/CONMIN  could  optimize  the  objective  function  satisfying 
the  inequality  constraints,  a  reliable  heat  balance  could  not 
be  obtained.  The  solution  to  this  problem  was  the  addition 
of  the  ALMM  option  to  COPES.  In  this  way,  the  method  of 
feasible  directions,  which  works  best  with  inequality  con¬ 
straints,  was  used  to  satisfy  the  inequality  constraints. 

The  multiplier  method,  which  works  best  with  equality  con¬ 
straints,  was  used  to  perform  the  heat  balance. 

The  results  of  the  case  studies  show  that  HEDSUP  will 
yield  reliable  designs  for  various  design  objectives  and 
problems  with  only  some  trial  and  error  application  of  the 
initial  Lagrange  multiplier.  However,  precautions  must  be 
taken  to  overcome  the  relative  minima  that  plague  this  design 


problem.  Table  1  shows  vividly  the  problem  of  relative 
minima  with  one  "optimum"  design  having  a  volume  over  100% 
greater  than  the  "true"  optimum. 

The  value  of  numerical  optimization  in  a  design  problem 
of  this  size  cannot  be  overemphasized.  For  example,  in  case 
study  1,  the  problem  is  taking  place  in  a  nine-dimensional 
design  space  and  intuition  on  how  an  optimized  design  "should" 
turn  out  is  quickly  lost.  Figure  20  helps  to  illustrate  this 
point.  When  only  varying  two  design  variables,  h  and  w. 

Figure  20  shows  that,  when  beginning  from  the  initial  starting 
point  with  a  Type  7  configuration,  a  design  satisfying  both 
the  equality  constraint  and  the  air  pressure  drop  constraint 
could  never  be  found.  However,  Table  2  indicates  how  the 
numerical  optimization  routine  has  varied  the  other  seven 
design  variables  in  order  to  shift  the  constraints  and  yield 
an  optimized  design. 

The  results  show  that  AHDOP  did  vary  the  surface  design 
variables:  D . ,  D  ,  i,  t,  s,  P*.  and  PT ,  in  order  to  produce 
an  optimum  heat  exchanger.  This  capability  of  surface  opti¬ 
mization  is  dependent  upon  the  use  of  the  Briggs-Young  and 
Robinson-Briggs  Correlations,  eqs .  (32)  and  (34a),  respec¬ 
tively.  The  reliability  of  the  correlations  as  compared  to 
the  "conventional"  method  is  questionable.  Actual  experimental 
data  for  a  particular  tube  and  pitch  will  always  be  the  most 
useful  in  predicting  pressure  drop  and  film  coefficients  of 
across  banks  of  finned  tubes.  However,  the  correlations 
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▼ 


1 


mentioned  above  cover  the  ranges  and  pitches  used  in  air¬ 
cooled  heat  exchangers,  and  should  therefore  be  sufficiently 


accurate  in  predicting  Hq  and  Apa> 


I 
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VI .  RECOMMENDATIONS 


In  addition  to  the  insight  that  this  investigation  has 
given  into  the  generation  of  an  automated  air-cooled  heat 
exchanger  design,  it  has  also  generated  an  awareness  of  this 
investigation's  shortcomings.  As  mentioned  in  the  review  of 
previous  work  in  this  area,  each  of  the  optimization  methods 
has  its  own  limitations;  none  is  completely  general.  Pre¬ 
sented  herein  are  recommendations  for  improving  upon  and 
furthering  development  of  HEDSUP. 

1.  HEDSUP  should  be  expanded  to  include  the  capability 
for  two-phase  tubeside  fluids.  Mott,  et  al.  [2]  discusses  a 
method  involving  two-phase  tubeside  fluids  that  would  be 
compatible  with  HEDSUP.  The  modular  design  of  ANALIZ  will 
aid  this  effort. 

2.  Research  in  the  area  of  numerical  optimization  using 
discrete  variables  would  benefit  HEDSUP  immensely.  With  dis¬ 
crete  variables,  the  design  of  the  exchanger  could  be  accom¬ 
plished  with  "off-the-shelf"  materials.  At  present,  the  use 
of  the  optimizer  is  restricted  to  continuous  variables  with 
continuous  first  derivatives.  The  ability  to  work  with  dis¬ 
crete  variables  would  also  eliminate  the  need  for  the  design 
matrix.  The  optimizer  could  optimize  for  type  configuration, 
fin  profile,  fin  material  and  tube  material. 

3 .  Additional  research  with  ALMM  is  needed  to  remove  the 
"artwork"  involved  with  choosing  an  initial  multiplier  and 
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thus  increase  reliability  and  hence  reduce  CPU  time.  The 
research  should  be  concentrated  in  two  areas:  1)  debugging 
the  new  optimizer  with  AFMM  and  2)  scaling  of  design  variables. 

4.  The  addition  of  cost  as  an  objective  function  would 
increase  the  attractiveness  of  HEDSUP.  Mott,  et  al.  [2], 
and  Fontein  and  Wass ink,  [18],  have  presented  much  useful 
information  in  this  regard. 

5.  Mechanical  constraints  such  as  tube  bursting  stress 
and  tube  vibrations  should  be  included  in  HEDSUP. 
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VII.  FIGURES 


Figure  1 
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NUMERICAL  OPTIMIZATION  TECHNIQUES 
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(JSADLE-F£ASIBLE  DIRECTION 
|/\,VF<X<!) 


s'!  FEACICLE 


Figure  4 


TUBE  PITCH/DIAMETER  RATIO 


Figure  5 


TUOE  PITCH/DI AhETER  RATIO 


0.0  0.2  0.4  0.6  0.8 


TUBE  OUTSIDE  DIAMETER  (INCH) 

Figure  6 
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CONTOUR  OF  CONSTANT  VOLUME 


TUBE  PITCH/DIAMETER  RATIO 


! 


Figure  8 
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TUBE  PITCH/DIAMETER  RATIO 
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Figure  9 
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Figure  10.  Cenf iguratien  ef  Air-Ceeled  Heat 
Exchanger 
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VISCOSITY  Cl.BM/FT-HR) 


VISCOSITY  OF  WRTER  (Rndrade's  &  Exp. 


VISCOSITY  (Ibm/ft 


CASE  STUDY  ONE 


CASE  STUDY  TWO 


OPTIMUM  DESIGN 


APPENDIX  A 


HEAT  EXCHANGER  DESIGN  USING  NUMERICAL  OPTIMIZATION  (HEDSUP) 

USERS  MANUAL 

In  order  to  execute  HEDSUP,  it  is  necessary  to  provide 
formatted  data  for  COPES  [25] ,  followed  by  formatted  data 
for  the  ANALIZ  portion  of  the  program.  This  section  defines 
the  data  which  must  be  supplied.  The  data  is  submitted  in 
segmented  blocks.  All  formats  are  alphanumeric  for  TITLE 
and  END  cards,  F10  for  real,  and  110  for  integer  data  when 
supplying  COPES.  For  ANALIZ,  formats  are  F14  for  real  and  110 
for  integer  data. 

While  the  COPES  user's  sheets  define  the  data  in  formatted 
fields  of  ten,  the  COPES  manual  does  provide  means  of  simpli¬ 
fying  this  task  through  unformatted  data  input. 

The  included  GLOBAL  CATALOG  defines  objective  functions, 
design  variables  and  constraints  along  with  their  location 
for  ease  in  compiling  the  necessary  COPES  data. 


GLOBAL  CATALOG 


GLOBAL 

LOCATION 

FORTRAN 

NAME 

DEFINITION 

i 

FMDOT 

Fluid  Mass  Flow  Rate,  lbm/hr 

2 

FLTMP1 

Inlet  Fluid  Temperature,  °F 

3 

FLTMP2 

Outlet  Fluid  Temperature,  °F 

4 

AMDOT 

Air  Mass  Flow  Rate,  lbm/hr 

5 

ARTMP1 

Inlet  Air  Temperature,  °F 

6 

ARTMP2 

Outlet  Air  Temperature,  °F 

'  7 

TUBEID 

Tube  Inside  Diamter,  in 

8 

TUBEOD 

Tube  Outside  Diameter,  in 

9 

FINHT 

Fin  Height,  in 

10 

FINTH 

Fin  Thickness,  in 

11 

FINSP 

Fin  Spacing,  in 

12 

PITCHL 

Longitudinal  Pitch,  in 

13 

PITCHN 

Transverse  Pitch,  in 

14 

BANKW 

Bank  Width,  in 

15 

BANKH 

Bank  Height,  in 

16 

VOLUME 

Volume,  ft 

17 

AREA 

Heat  Transfer  Area,  ft^ 

18 

FHP 

Air  Horsepower,  HP 

19 

PPA 

Airside  Pressure  Drop,  psi 

20 

PPW 

Tubeside  Pressure  Drop,  psi 

21 

DELQ1 

Q  -  Q5,  BTU/hr 

22 

DELQ2 

Q5  -  Q4,  BTU/hr 

23 

DELQ3 

Q4  -  Q3,  BTU/hr 

24 

DRATIO 

Df/°o 

25 

TUBTH 

Tube  Wall  Thickness,  in 

26 

TOUCHN 

Tip-to-Tip  Clearance,  Transverse 

27 

TOUCHL 

Tip-to-Tip'  Clearance,  Longitudinal 
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GL03AL  CATALOG  (Cent.) 


28 

PROFH 

29 

PROFC 

mmmammmmmmi 

30 

ARG5 

1  -  1 2Apo/pooJ  .  . 

31 

QDOT1 

Q*.  BTU/hr 

32 

DELPHI 

Q-Qmax*  ^dians 

33 

QRATIO 

Qs/Q 

34  ' 

ARG7 

V2pr. _ 

35 

ARG8 

Df/2PL 

36 

VROWR 

Number  of  Tubes  per  Vertical  Row 

37 

DELSFF 

Negative  of  Free  Face  Area 

38 

THETA 

arcsin  (ARG7) 

j 
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DATA  BLOCK 


DATA  BLOCK 


Print  above  plus  each  proposed  design  vector ,  objective  function  and 
constraint  values  during  the  one-dimensional  search. 


FIELD 


AD-A096  350  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA 
HEAT  EXCHAN6ER  OPTIMIZATION. (U) 

SEP  80  C  P  HEODERICH 

UNCLASSIFIED  _ 


2  of  2 


FIELD  CONTENTS 


of  the  objective. 


DATA  BLOCK  F  OMIT  IF  Nl)V  -  0  IN  BLOCK  B 


DATA  BLOCK  It  OMIT  If  MDV  =  0  IN  BLOCK  B 


DATA  BLOCK  I  OMIT  IK  NOV  =  0  IN  BLOCK  B.  OR  NCONS 
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DATA  BLOCK 


ANALIZ  DATA 


3)  Comment  cards  are  no  longer  permitted. 


DESCRIPTION:  Airside  Parameters 


DESCRIPTION:  Tube  Geometry 


DESCRIPTION:  Tube  Arrangement 


DESCRIPTION:  Integer  Parameters 


DESCRIPTION:  Miscellaneous  Parameters 


DESCRIPTION:  Penalty  Parameters 


DESCRIPTION:  Optimization  Parameters 


DESCRIPTION:  Augmented  Lagrangian  Function  Method  Data 


COPES  DATA 


DATA  BLOCK  F  -  COOT 


DATA  111.0CK  G  -  OMIT  IF  NDV 


SCAL1  BU  SCAL2 


DATA  BLOCK  I  -  CONT 


DATA  BLOCK  I  -  CONT 


DATA  UI.OCK  V  -  COl'LS  END  OF  DATA  CARD 


APPENDIX  C 


SAMPLE  OUTPUT  FROM  COPES 


aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa aa aa aaaa aaaaaaaaaaaaaaaaaaaaaaaaaa aaaaaaaa aaaa 


03/16/80  13. 5*. 31 


FILES  SAV 

1 


DATA 


T  1 


NAVAL  POSTGRAOUATE  SChOOL 


^cccccc  • 

£ 

c 

c 

ccccccc 


00 COCO Q 
0  0 

8  8 

U  0 

0  0 

occcooa 


ppppppp 
p  p 

'ppppppp 

p 

p 

p 


|ESEEEE 

^EE: 

E 

E 

EEEEEEE 


|SSSSSS 

Issssss 

s 

s 

sssssss 


CONTROL  PROGRAM 
FOR 

ENGINEERING  SYNTHESIS 


1 

0 


title 

AIR  HEATER  (DESIGN  -  OPTIMIZE  VOLUME. 
CARC  IMAGES  CF  CONTROL  JA  TA 


CARO 

IMAGE 

1  ) 

AIR  HEATER  DESIGN  -  OPTIMIZE 

2) 

2.9 

2) 

2  9 

3) 

i  BLOCK  C 

4) 

0.90.0.10 

41 

0  30 

0 

51 

i  9 LOCK  0 

6) 

0. 

6) 

0. 

71 

0. 

i 

1 

0. 

I 

1 

$  BLOCK  E 

< 

1 

9. 1 6 »- l . 

« 

1 

9  16 

-l. 

1C 

» 

t  BLOCK  F 

11 

l ) 

.232,2. J25,. 6736 

I 

f 

.  >232  ,  2.325 

.6736 

1< 

!  1 

.25, 2. 5.. 7 201 

1. 

I) 

.25  2.5 

.7201 

10 


St  V 

m 
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DATA  11 

.0625, 1.0*20,. 16 
.Oo?a  1  . 

0+2  0 

.  16 

14  ) 
14  1 

.01, .021/. .0235 
•  0  1  . 

0217 

.0235 

l* ) 
IS ) 

.03,. 126, .03 
.08 

.125 

.08 

16)* 

16) 

0. , 9. ,  1.04 

0. 

4. 

1.34 

IT) 

17) 

0.  ,  4.  ,  1.04 

0. 

4. 

1.34 

18) 

18) 

0 . • 500 . ,236. 

0. 

500. 

29  5. 

11) 

0.  , 500. ,342  . 

342. 

19) 

0. 

500. 

20) 
21  I 
21) 

S  BLOCK  G 

U' 

7 

22) 

22) 

2, a 

2 

6 

23) 

3,9 

23) 

3 

9 

24) 

24) 

4,10 

4 

10 

25) 

5,11 

11 

25) 

5 

26) 

26) 

6,12 

6 

12 

27) 

271 

7,13 

7 

13 

28) 

23) 

8,14 

8 

14 

29) 

29) 

9,15 

9 

15 

30) 

31) 
31) 

ill 

t  BLOCK  h 

u 

11 

33) 

l.,0.,2.5 

33) 

1  • 

0. 

2.5 

34) 

34) 

35) 
35) 

25  25 

.  01  9,0., .18 
.0  13 

0. 

.18 

36) 

36) 

20,27 

26 

27 

37) 

37) 

-1.0+20,0.,  0. 
-1.0+20 

0. 

0. 

39) 

33) 

391 

39) 

19 

19 

-1.0+20,0.  ,  .  0722 
-1.0+20 

0. 

.0722 

40) 

40) 

41 ) 

20 

20 

—  1 . 0+  20,0.,. 14 

N  AVAL  f*3SI  \XAOUATE  SCHOOL 
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FILE:  SAV 

DATA  T1 

41! 

-l  .0*20 

0. 

.14 

42) 

30 

42) 

30 

43  1 

•  7 , 0  * ,  1.0*20 

43) 

.  7 

0. 

1.0*20 

44) 

38,38 

44)* 

38 

38 

45  ) 

O.OtJ.0,1. 3  ,0. 

45) 

0.0 

0.0 

1.3 

•*6  ) 

34,33 

46) 

34 

35 

47) 

■1,0  ,0.  ,  1. 

47) 

-1.  J 

◦  . 

l  . 

4d ) 

36 

48) 

36 

49) 

2.  ,0.  ,  1.0*7 

49  ) 

2. 

0. 

1.0*  7 

50) 

37 

50) 

37 

51  ) 

-l • 0*2  J ,0 . ,0. 

51) 

-1 . C*20 

0. 

0. 

52) 

33,33  ,2 

52) 

33 

33 

2 

53) 

1.  ,  .  1 , 1.0*20 

53) 

1. 

.  1 

i  .0*20 

54! 

ENO 

1  TITLE 

• 

AIK 

HEATER  DESIGN  - 

OPT  I. MIZE 

VOLUME 

NAVAL  PO  STijR  AGUAT  E  SCHOOL 


0. 


2 

9 

J 

0 


CALCULATICN  CONTROL,  NCALC 
VALUE  '16*0  1  40 

1  SINGLE  ANALYSIS 

2  OPT  I'lIZAl  1  u  J 

3  SENSITIVITY 

A  TwJ-VAk U JLZ  EJECTION  SPACE 

5  0PTIM.J1  itNSITl/ITY 

6  APPROXIMATE  OPT IH It  AT  1  ON 


CCNTKCL  PARA1ETERS, 

C ALCJLAT I C  V  CCNTPOL,  NCALC  = 

NUMBER  JF  ,|  jJU  JtjlGN  VARIABLES,  NOV  * 
INPUT  I  Nr  HIM  AT  ION  PRINT  COu=,  IPNPUT  = 

OEBUG  PRINT  CODE,  IPODS  * 


*  *  OPTIMIZATION  INFORMATION 


global  VARIABLE  ivUMJER  OF  OBJECTIVE 
MULTIPLIED  (NEGATIVE  INDICATES  MINIMIZATION! 


*  -0 


16 

.IOOCE  01 
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FILE 


s.w 

DA  T4 

T  1 

NAVAL 

p'JSTgraouatc 

SCHOOL 

12 

34 

35 

0 

-0.  10  0  00  E 

01 

o.iooooe 

0.  20  0  )0H 

01 

0.100 00 5 

01 

16  * 

3b 

0 

0 

0.  2  )00 JE 

)l 

01 

O.lOtJJt 

Of 

IB 

J  7 

0 

0 

-0.  1  no;)E 

16 

0.  1  l  U  J  J  C 

1.4 

0  .0 

19 

33 

Jo 

2 

0.  IOOOOl 

01 

0 . 1 0  OOOt 

JO 

0. 1  10 JOE 

u 

TCTAL 

NUMBER  OF 

CUNiTRAlNEO 

PA RAM* IE  RS  * 

13 

*  •  ESTIMATED  DAJA  STLRAGE  REDO 1RCMENTS 

REAL  INTEGER 

INPUT  EXECUTION  AVAILABLE  INPUT  EXECUTION  AVAILABLE 

106  550  SCOJ  7f>  146  1000 

A I  R-CCQLED  MEAT  EXCHANCER  0>T  l  MUAT  I  C.M 


INPUT  DATA 


TUBESIDE  PARAMETERS 

MASS  F L Q'.<  RATE*  133333.3125  L3M/HR 
INLET  TcMPEKATJ.it*  200.0000  3E3  P 

OUTLET  TE  M°c  RAT  UK  E*  125.0UJJ  OSG  F 

SPECIFIC  HEAT*  i.ojoo  BIU/L  dM-F 


AlRSlOE  PARAMETERS 


MASS  FLCW  RATE* 
inlet  t empce aturc 
OiJTLeT  TEMPcRATUr 
SPECIFIC  HfAl* 
INLET  PRESSURE* 


1190476.0000  IBM /HR 

95.  JO  00  >EG  F 
*  1 30.0l00  JEG  r 

0.24 GO  bTU/LrtM-F 
14.0000  PSl 


TUBE  GEOMETRY 

TUBE  iriilOE  0 1  A .  *  0.8700  INCHES 

TUbt  CUTS  I  j t  JIA.*  1.080)  INCHES 

FIN  HEIGHT*  0.4600  INCHES 

FIN  IhICiLNESS*  0.0190  INCHES 

FIN  SPACING,  CEITER-T J-CtNI EH*  0. LllO  INCHES 


TUBE  ARRANGEMENT 


130 


FI IE l  SAV  OAT A  Tl  NAVAL  POSTGRADUATE  SCHOOL 


TRANSVERSE  PITCH*  2.1250  fNCHFS 

LONGIIUj  I  NAc  PI  I.ll*  2.1253  INCHES 

BANK  HEIGH*  46.0000  INCHES 

BANK  RUTH*  238.0000  INCHES 


INTEGER  PARAMETERS 

TYPE  CF  CROSS  FCUv,  ARRANGEMENT*  7 

TYPE  CF  FIN  p  \HT  I  L  c*  1 

NUMBER  OF  Rj.iS*  4 

NUM4E  R  OF  PASSES*  2 


MISCELLANEOUS  VARIABLES 

GIVEN  FEAT  TRANSFER  RATE*  1000 3000 . 3 003  8TU/HR 

THERMAL  CONDUCT  W  i  U  OF  TU3  =  MATERIAL*  220.  0000  3TJ/HR-FT-F 

THER/vAL  CO.VJUCTl  Vi  TY  CF  FIN  MATERIAL*  U8.0000  cJTU/HR-FT-F 

GAS  CONSTANT*  53.3400  FT-LBF/LHM-R 


penalty  parameters 

12. CCCOCOOOO  0.0  0.0 


OPTIMIZATION  PARAMETERS 


0.0 


0.0 


0.0 


0.0 


*  CC  *  0. 100 JOE  02  C MOL T  *  0.20000c  01  CCMAX  *  0.1030DE  04 
KOUNT*  l 


•  •  COMIN  0EIEC1S  INITIAL  Xin.GT.VUBU) 

Xllt  *  0.23  >jc  -  0 1  V.Jim  *  0.2  173E-J  l 

xm  is  set  c-joal  to  vjji  ii  for  i  =  4 

KCUNT*  1  JEL*  0.  7  74.4  4F  01  CBJ*  3.753  10E  02  04  01  *  J.  723  C9E  02 

—0. 435  lOt  00-0.^25'ioi:  JO-O.  2>)S  S-u  3-  3  .  9000  OC  oO-0  . 23553c  01 
— 0. 8  44  04  E  —  0 1  -  ') •  jO  '3  2  >  (. -  J  i  —  3  •  f>9  79c  _  —  01  —\j  •  4  14  '  4  E  3J-3.o74J4E  01 
— 0.47/C1F  JO-,1.  1  32.- 7E  J  l  -  3.  .>7  i  2  *  F  JJ-0.14932F  ol  -0 . 5  04  7 J  E  00 
— 0.  1 0044  E  03-0.  ■>  yf;03t  00-0.  1  6  7  73  E  06  0.774446  00 


0.0 

0.  0 

0.3 

0.  0 

3.0 

0.0 

0.0 

0.  0 

0.0 

0  .0 

0.0 

a.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.  0 

3.77434E  01 
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FILES  SAV  OATA  T1  NAVAL  f»3  ST  UK  AOUATE  SCHOOL 

KOUNT *  2 

KOUNT -  2  DEL*  0.22126F  01  UiU*  0.71  452C  02  0001  *  0.72922E  02 

-0.454046  JO— J  .  4  IM  3*>  t  JO-'). 132  J  46-01-0.  3444)1 1=  OJ-J. 191205  01 
“3  •  17flr>  lf-0'j-  J.  3  7  72  1L - 02-0. 424 2 l E- 01 - > .4 1  3  4  36  03-J.OJ24Q6  1  l 
-0. 447/Oh  00-J.  16 > /vt  Jl-0.3  <2416  CO-J.  150UG  JI-0.5JJ006  JO 
-O.lliSit  Oi-0.‘>)tsoL  JO-O.  lJ-»0  3fc  Oo-O. 221266  00 


0.0 

0.0 

0.0 

J.O 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0. 55358E 

01 

KOUNT3 

3 

KOUNT  3  3  DEL3  C.7C331S  00  00J=  0.71S79E  02  OBJ  1  3  0.722436  02 

-0.454776  00  -0 . 4  l  <3  )4  t  00-0 . 3  36626  -  01-0.  0^66  3  E  00-0.106206  01 
-0.  260  04  6-0  l -0.  1  ‘.>0  73c -  01  -  0.J71  70  6- 02-3  .4141  LE  JO- J  .on  loc  )  1 
— 0 .  4  Hi  1  ■> c  J 0—0 .  1  J -» 1  6  J 1 -  6 • 3  vuo  / 6  00— J •  14° 06c  01—  O  •  5  J  l  o  5  E  JO 
-0.11315E  03-0. 9 J94d fc  JO-0. 164986  06-0.700  1  lc-01 


0.0 

0.0 

0.  0 

3.3 

0.0 

0.0 

0.  0 

0  .0 

0.  0 

0.0 

0.0 

0.0 

0.0 

0.0 

0  .0 

0.0 

KOUNT* 

0.0 

4 

0.0 

0. 48320  E  01 

KOUNT=  4  DEL3  0.1C33OE  01  03J=  0.71241E  02  J0J1  3  0.71464E  02 

-0.454456  00-0. 41820c  00-0. I rf  )  O^c- Cl- 3 . 9  98 1 -6  00-0.13535E  01 
-0.8  70236-  3 4-0 .321  306-02-0.  0  4  0*.  5 E- 02-3  .  •*!  3  >2  S  00-0. 06  /  3  36  01 
—  0.4  94  3o  c  uO-O.lollUl:  JI-0.33J49E  C0-J.150J0E  J1-0.60001E  JO 
-0.11JC9E  J3-J.S!4sMt  JO— 0.  16  3  96  E  06-0 . 51  65  IE- Jl 
0.0  0.0  0.0  0.0  0.0 

0.0  0.0  0.0  0.0  J.O 

0.0  0.0  0.0  0.0  0.0 

0.0  0.0  0.0  0. 37990  E  01 

KOUNT 3  5 

KOUNT=  5  OEL=  0.100V3E  01  03J=»  0  .71294E  02  J3JI  3  0.  7  14596  02 

-0.45-.51E  OJ-O. 413206  00-0  .  I  6-  S4E-  Cl  -  0. 3  9R3  5:  00-0.  135J3E  Jl 
-0.566  94  6-0  i-0 . 3 3 o 02  6-02 -0. 3 49 7  l  S- 02- J . 1 1 342 E  00- 0.63 7  3  36  0  1 
-0.494366  00-0.  IM  IOC  Jt-0.3rf.J09t  00-0.  I  50  Joe  J1-0.5JJ03E  JO 
-0.11 30cc  0 3  - J . 94493  fc  00-0.  16  39  /E  06-0 . 5024 CE-Ol 


0.0 

J  .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0  .0 

0.0 

3.0 

0.0 

0.0 

0.0 

0.0 

0. 27942  E  01 

KOUNT* 

6 

KOUNT*  6  OEL=  0.94123E  00  C3J 3  0.  71  3  38  6  32  03JI  3  0.714486  02 
-0 .454656  00—0.4 131  3  6  30-0.  1 2 7 232- 3 1- J . 3 59 73 E  03-9.133996  31 
-0.1  70  17  fc-3  2-0 . 3  70  3  *6- 02— 3.08  J<..;-02  -  )  .4  1  3  i  JE  J3-0.oj733c  01 
-0.49436E  0.3- J .  1  o  l  1  J  6  J  l-  3.  3<-9  OJE  Co-0. 14499?  Ol-J.pJOlJE  30 


-0.1  I309t 
0.0 

03- 0. 449406 
0.0 

Jj-O.  16  399E 
0.0 

3  0-  J  .  4  730-.E-  01 
0.0 

3  .0 

0.0 

0.0 

0.0 

0.0 

J  .0 

0.0 

0.0 

0.0 

0.0 

3  .0 

0.0 

0.0 

0.0 

0.18529E  01 
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FI IE :  SAV 
KOUNT  * 


7 


DATA 


U 


NAVAL  POST  GRADUATE  SCHOOL 


KUIM*  7  DEL*  0.103196  Ot  URJ*  0.714026  02  T!3Jl  *  0.71445E  02 
-0.454666  UO-J.'t  l  9  18C  JU-J.t  1355E-U-0.  Js.l-sr>6  JO-O. 1(56036  01 
-0  .2  12  3  7L-  )2  -  J.  3  J  /■.*  4l-J2-  )  .  39  1 2oE-  02-  0  .•*  1  10  V:  OJ-).oj/ JJE  3  1 
-0.494  V>t  JO-O.lLiUjL  Jl  -  0.  3  a  >00 1  JO-J  .  149  74  E  Jl-J  .5U0lJt  00 
-0.1  1  3  0 9 L  03- J.  >  >994  t  JO-  J  .  I  6‘»  OUfc  06-0. 457956-01 


0.0 

0.3 

0.0 

J.3 

3.3 

0.0 

3.0 

3.3 

0.0 

3  .3 

0.0 

0.  3 

3.  J 

0.3 

0.0 

0.0 

KOUNT* 

0  .J 

3 

3.0 

0.213666-01 

KOUNT *  9  OEL*  0.  17517E  01  04J  =  0.71433E  02  JBJl  *  0.71445E  02 

-0.45452L  OJ-0.4LH10L  JO-O.l  13106-01- ).  H9.1B7E  OJ-O.H606E  01 
— 0.2 435 4E—C2— 0.41  4  26  E-J  2-0  .  J-  1  Z  IE-  32-0.41  j  >4E  00 -0  .a  5  7  3  3E  .01 
— 0. 494  3o  E  00-0  .  l--,l  lot  JL  -0.3»4  JOE  00-3  .  1 4  >49  E  Jl-J. 50016c  00 
— 3.1UIJL  wJ-J.>>4>Jt  JO-O  .  1  6 •»  0J 6  06-'J  .  4  3  7  >  J£— 0  l 


0.0 

3.0 

0.0 

0.0 

0.0 

0.0 

3.3 

0.0 

0.  0 

0.3 

0.0 

0.0 

KOUNT* 

0.0 

3.3 

9 

0.0 

0.0 

0.0 

-0.  173066  31 

0.3 

KOUNT*  9  36 L*  C.16433E.01  JaJ =  0.715556  32  U8JI  =  0.71446E  02 

—0.45449c  Jj-0.41o20E  JO -0 . 1 12  &o£- 2 1  -  0.  d  43'J  7  6  OJ-J.ldbJaE  31 
— 0 .264596-02-0.  4  3  344  fc- J  2-3 . 3a  6  05C- 02-0  •  1 3  >4£  00-3  .66  1  3  3E  31 
—0.4  94  34  fc  JJ-J.l6lljt-  jl  -0.3  6  3  00  E  00-3.  144046  Ji-J.iJOlbE  00 
-0. 11310c  OJ-0.  9 -*9  40  6  JO-0 .1  6->C5E  Jo- J . 424586-0 l 
0.0  0.0  0.3  0.3  0.0 

0.0  0.0  3.0  0.0  0.0 

0.0  0.0  0.3  0.0  0.0 

0.0  0.0  3.0  -0.342396  31 

KOUNT*  13 

KOUNT*  10.  DEL*  0.25459E  01  U3J*  0.71530E  02  OBJl  *  0.71599E  02 

-0.45272E  03-J.4l.U16  JO-  3.  1  3624c -01-3.  64965  E  Jj-0.l34d46  31 
-0. 14077C-Jl-0.l630lL-3l-0.55 76-6-02-3.414116  uJ-3.a4702E  31 
-0.  49460  6  JO-J.  lilo-.t  01-0.36  >24E  33 -3  .  144  )  2  E  0  1  -3 . 6  JO  e4E  00 
— 0. 1 1  J4'.E  03-0. 4V4>Jc  J  J- J  .  lo5  22E  J6  0.  318246-31 


0.0 

0.0 

0.0 

J.O 

0.0 

0.0 

J.O 

3.0 

0.  0 

0  .0 

0.0 

3.0 

0.  J 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.883326  00 

KOUNT* 

11 

KOUNT*  ll  3EL*  1.25754E  01  OHJ*  0.715396  02  33JI  *  0.715966  02 
-0.45271c  JJ-O. 413916  )0  -  0  .  I  31  J  56  -  Ol -3.  ei  455  75  J0-J.15664E  01 
-0. 141495- J 1-0. 152 /5c -31-3.54 1996-32-3 .4141 l 6  3 J- 3.65 701 E  31 
-0.494.6lE  JJ-0 •  lo  1 0  1  c  J  l  -0.3  4  >2  4  c  0J-3  .  1494  2E  3  1-0.6  30(146  JO 
-3.1 1 J44£  C3- J . 9  >4  >o  L  30-3.  lo522c  06  3.321926-01 


0.0 

0  .0 

0.0 

0.0 

0  .3 

0.0 

0.0 

3.3 

0.0 

0.0 

0.0 

J.O 

0.3 

0.0 

3.0 

0.0 

0.0 

3.0 

0.169246  01 
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FILE:  SAV  OATA  Tl  NAVAL  POSTGRADUATE  SCHOOL 

KOUNT«  12 

KGUNT*  12  1)6 L  =  0.25936b  01  IJfiJ  =  0.71691E  02  JAJl  *  0.71594E  02 

“0  .9  52  711  JO- J.4l  9521  JJ-O.  12JO/L- JI-0.39H/11;  JU-J  .  L  34  I J  L  Jl 
-J. 141 S4C-J  1-0.  16  2-j  7  t  -  H— ).  an  J  *  3E-02-9.4  1411  c  )J-J  .6)  7006  01 
-0.49461t  JJ-<).  1-jIOM:  Jl-O.JJ  92:>6  00-0  .  I  -.492  6  01  -J  .  'J  JO  »  5  E  JJ 
— J.1134j6  J  j- J  .  9  94  9  J  6  JO-J . 1  to 22 C  06  J.3242JE-JI 


0.0 

J.  0 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

3.0  0.0 

0.0 

0.0 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

0.42860E  01 

KOUN T* 

13 

KCU  NT  *  13  OEI.  =  3.17131E  01  OBJ*  0.71739E  32  30J1  *  0.71  776E  C2 

-0.451766  JO-J  .  ♦  1  92  96  JJ-0.20nj4  E-Ol-J  .09  7  )  4  E  00-0.135066  01 
-0. 177J3C-0  1-0.  U3J  .U- Jl—J.l  14  32  2-Cl- J.41-.  )  7t  JJ-J.  t>  5o  7  1  £  H 
—0  •  4  9  •«  4  -» E  )J-0.1*>1  Jar.  01  -0.36  J*»‘J  E  JJ-0.149396  Jl-0.5J10oE  JO 
-0 .113  726  JJ-O.  /  jvSoc  JJ-0.l6-.47f:  C6- 0.1  0733  c— 01 


0.0 

0.  0 

0  .0 

3.0 

0.3 

0.0 

0.0 

0.  0 

3.0 

0  .0 

0.0 

0.0 

0  .0 

0.  0 

0.0 

0.0 

0.0 

0.0 

3.25679E 

01 

KOUN  T= 

14 

KCUNT  =  14  DEL=  -3.1c953E  01  00  J  =  0.71  757c  02  08J1  =  0.  71  775E  C2 

—0.451786  J  J— 0.419236  ui)-  0. 2  04  6  3  E  —  0 1  -3  . 3  9  7  9  6  E  00  —  J.IjoJoE  01 
— 0  .  I  7  7  4Cc-  J  l  -0  .  U50.C-  Jl-O.  1 1  5  02C-  01-).  4  1-0  7E  JJ-J.o567lc  Jt 


-0. 49484E 

0  J  —  J  .  1  o  1  0  j  c 

Jl-O. 33 >-4fe 

JO-J . 1 49S9E 

01  -  0  .  a  0 1  Jo  c  JO 

— 0. 1 13  72  t 

0  3-0.  9 9  9  9  3c 

00-  0.1  lit  97E 

J6- J.  1  0595 E- 

01 

0.0 

0.  J 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3.0  ■ 

0.0 

0.0 

0.0 

0.0 

Q.  0 

0.0 

0.0 

0.0 

0.0 

3 . 8  726  CE 

00 

KOUNT* 

15 

KCJ  NT  *  15-  0EL=  J.  16739E  01  OOJ=  0.71774E  02  JBJL  =0.7V774E  C2 

—0.451 7b  6  0 J-0 .  4 1 42  9C  00-0. 2  Cl  1 7 E- 0 1- 3 . o9 7 ) 7 E  JJ-0.1363oE  01 
-0.1 78  35  t-0 l - J .  13  50.C- J 1-J. 1 1  5  62E-01-J.414J7:  00-0. o -5  7  12  Jl 
— 0.  4  44tWE  00-0 . 1  -j!  J6  c  31— J. 39949c  OJ-J .149396  01-J.501J7E  00 
-0.U372E  J  3—0 . 99  4  4  J  t  J0-J.16497E  Jb-O.  10*.6  2E-01 


0.0 

J.O 

0.  J 

•3.0 

0.0 

0.0 

0.0 

0.0 

3.3 

0.0 

0.0 

0.  J 

0.3 

3.  0 

0.0 

0.0 

0  .0 

0.0 

-3. 80129E 

00 

KOUN  T* 

16 

KCU  NT* 

16  JEL  =  0 

I.3J295E  01 

QOJ=  0.71799E  32  3BJI 

=  3.  717736 

-0.451  7JE 

00-  J.  41  92  9  6 

J0-0.2u2?6E- 

H- J  .3  9  79  3  F 

JO  -u  •  l  ^cidoE 

Jl 

-0. 1 73O0E 

-Jl-J.lJ-.93i.- 

■J  i-J  .  1  16  tVE- 

0  L  —  U  •  4  l  *  )  7£ 

Jj- 3.566  7  16 

91 

-0. 494  J4£ 

30-J.loi  JJt 

J1-J.3J9  4.C 

JJ-3 . 1 4J596 

0 1  —  0  •  *3  J  i  0  7  fc 

JO 

— 0 .11371b 

J  3  - J . 9 9445  L 

JO- J  •  1  0-t  S  76 

20- J.  1  04  35  E- 

■ji 

0.0 

0.  J 

0.  J 

0 . 0 

3.0 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0 

J.O 

J.  0 

0.0 

3.3 

0.0 

0  .J 

0.0 

-0. 4  130  BE 

01 
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FILE:  SAV  DATA  II  NAVAL  PO  ST3*  AOCIATE  SCHOOL 

KQ  UN  T  =  1/ 

KOUNT*  17  OEL  =  0.21240E  01  D9J  =  0.71S03E  02  U3JI  =  0.71924F  02 
-0.46134C  JJ-O.  419-WC  jO- J.  2221  JL-Jl-J.  d4/7;»E  0 J-0  .  i Jo  y2  E  11 
-0.  l47d2t-0  l  -  )  .  Ibbl  3l-JI-0.  1  15  2  1C-  Jl- J.  4  1  '.l  ?C  0O- J.6b 340c  01 
-0.4T.W2F  00- ).  l<»  U-»l  U-J  .  J.i  -U7F  CJ-J.  1  49s;  d  I:  Ol-J.aJllSE  00 
—0 .113311  0 3  —  0 •  i'i'i  f-i L  JO- 0.  lob  2b c  Os  0.66  3  7  oE  —  02 


0.0 

0.) 

0.0 

0.0 

J  .0 

0  .0 

0.0 

0.0 

0.  0 

0.0 

0.0 

0.0 

0.0 

0.0 

0  .0 

0.0 

0.0 

0.0 

-0.200a  7E  01 

KOUNT-*  Id 

kount* 

18  ne  L= 

.0.21J53E  01 

OBJ*  0.71817E 

02  0UJ1  *  0.71823F  02 

-0.46134E  JJ-O. 419. 72  JO-0 . 22 1 32c- Cl -0. b 9 77 9c  OJ-J.Wo'WE  01 
-0. 199 J  1C- J  1-0 . I yol  3  l-JI-O. 1 15  l  ic-Jl-O. 4 14 12 C  00-0 .05  64 Os  01 
-0.49442  E  0')-).Ul.)*l  J  1  -0 . 3  *J  l‘j  7l:  CO-O  .  144  H  F  0 1-0.4011 oE  JO 
—0. 11331c  O-J-J.VVbVjL  00-0.  lo5 20 E  06  0.667456— 02 


0.0 

d.i) 

J.O 

J.O  J.O 

0.0 

0.0 

o.o 

0.0  J.O 

0.0 

0.0 

0.0 

0.0  J.O 

0.0 

0.0 

0.0 

0.12912=  JO 

KQ  UN  T  3 

19 

KOUNT=  10  OEL  =  0.42O76P  01  03J  =  0.71837E  02  0301  *  0.71822E  02 

-0.45134E  0J-J.410-.6l  00-0.  2  20  706-  01 -0.  8 977  96  JJ-J.Ui.v2E  )1 
-0.1  13  lOb-0  1-J.  16b 02 6- Jl-O. 115 l 5  6- Jl-O. 41412C  0J-u.55oe.0c  OL 
“0.4  9b 02 E  oO-O.lolJ-.E  J l -0 . 3 89 5 7E  CJ-J. 145886  J1-0.5JU3E  00 


-0. 1 13d  l  c  03 
0.0 

- J. S599jF 

0.0 

JO— J.  165  24  E 
0.0 

06  0. 67151 E- 02 
0.0 

0  .0 

0  .0 

0.  0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

KQ  UN  T  *  20 

0.0 

0.0 

0 .44268  E  01 

KOUNT*  20  OEl=  0.4J011E  <)1  OQJ=  0.71  366c  02  U3J1  =  0.71821E  02 
-0.461346  00- J . 4 1 9bo E  JG-0 . 220 JeC- J1 -0. 3 97du =  JO-O.U692E  H 
-0.  1  Jd  03  E-0  1-0 . 145  7J  c- Jl-O.  1 13  l?c-  Jl-o.  -.141  2  E  J0-0.o3oa  IE  0  l 
—  0.4  94  92  E  JO— J  •  16  l  J  t  c  )  1- J  .  3  'i9  b  76  GO-O.  I  493  4  E  0 1  -  J  .b  J  1 1  3  E  00 
-0 . 1 13d 1 E  0J-0.S9S9JI  JO-O.  16329c  Oo  0.672U4E-02 


0.0 

0.0 

0.0 

0.0  J.O 

0  .0 

0.  0 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

0 . 87273 E  Ol 

KOUNT  = 

20  061=  0. 

4  30  HE  01  03  J 

=  0.  71994c  02  3BJ1 

OPT  I.-l  12  AT  101  dtSJLTS 

OBJECT  WE  FUNCTION 

GL06AL 

LOCATION  15 

FJ  NCTJ  cn 

VALUE  0.71321C  02 

DESIGN  VAd TABLES 


FILE  :  SAV 

10 

1 

2 

3 

4 
6 

* 

8 

9 


DATA 

0.  V. 

NC. 

1 

2 

3 

4 

5 

6 

7 

8 
9 


T  l 

GLOBAL 
VAR.  NU. 
7 
5 
9 
10 
U 
12 

13 

14 

15 


NAVAL  POSTGRADUATE  SCHOOL 


LTUE  R 
•lU'iO 
0.2320 jE 
0.  250  JOE 
0.  n250.)C- 
0.949946- 
0.  799J9E- 
O.C 
0.  0 
0.0 
0.0 


Ov) 

00 

•01 

•02 

01 


V  AL'JE 

0. 5)H))€  JO 
0.  57554E  JO 
0.  1 2  9d  )b  OC 
O.  2)7  I  >E-Ol 
U. 79 04 )E -01 
O.dJlicE  JO 
0.  1022  1c  01 
0.  136075  0  3 
0.2360  )t  03 


01 

01 

lu 


UPPER 

nou.\c 
.2  )2 8  C! 

.  25)005 

.  1 100  V'- 

.21  7)J_-')l 

.12500b  JO 
.39908=  01 
.-.OOJOc  01 
.<.9895=  03 
.494946  03 


DESIGN  CONSTRAINTS 


ID 


GUJdAL 
VAR .  NO. 


LOWER 

aUJNO 


VALUE 


UPPER 

BOUND 


1 

24 

0.  1 00  COE 

01 

3.  L  45  L  36 

01 

0.2  53C0.5 

01 

3 

25 

0.  ldOJJE-Ol 

0.1R"W7E- 

01 

0*1$ J  OOE 

00 

5 

26 

-0. 1 100  jE 

16 

-0. 1 rtCM  2 1 

JO 

0.0 

6 

27 

-o.  uuooe 

16 

-  0 . 1  9d  ).■}£- 

J2 

0.0 

7 

19 

”0.11 00) E 

L  6 

3.7  10  J  35- 

jL 

0. 722CCE 

-01 

8 

20 

-0.110(136 

16 

0  .  t  2.1 5  9c 

0) 

0.  1 40 CUE 

00 

9 

30 

o. Tooyjt 

00 

0 .909RRE 

00 

o.  i  Door 

15 

10 

oa 

0.0 

0.  o566  JE 

JO 

0. 130)  JE 

31 

12 

34 

-0.  1)0 00  £ 

0  1 

0 . 61043= 

00 

0. IOOOOE 

3  1 

14 

35 

-0.10J03E 

01 

0 . 4  *30  2  c 

JO 

o.iooor; 

Jl 

16 

36 

0.2J0  0JE 

01 

0.  2  2  96 jE 

03 

0.  1  JO JOE 

3  3 

13 

37 

-0.1 1000E 

10 

-3  .  1  652  95 

05 

0.0 

19 

33 

0.  IOOOOE 

01 

0.  99933  £ 

JJ 

0 . 1  10 JOE 

16 

OUTPUT  DATA 


TUBESIDE  PARAMETERS 
NASS  FLOU  RAT==  133333.3125  IBM/HR 

inlet  temperature*  200.0020  oeg  f 

OUTLET  TEMPERATURE*  12  5.0  3  30  JES  F 

SPECIFIC  FEAT-  1.0000  QTU/IB'1-F 


AIRSIDE  PARAMETERS 
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FILES  SAV 


DATA 


n 


NAVAL  POSTGRADUATE  SCHOOL 


MAS  S  FLOW  HAT  P»  U?047S,.00J3  L'lM/HR 

INLET  IEHs»Ek  ATiWs*  U5.00CJ  l)tj  F 

OUTLET  I E  f*i*L  H.%  TtjHli*  130.0  5  10  JLG  F 

SPECIFIC  HE  A  I  *  ■  0.2400  iJTU/lU'l-P 

INLET  PRESSURE*  U.JOOO  PSl 


TUBE  GEOMETRY 

TUBE  INSTDF  DIA.*  0.5180  INCHES 

Tuae  uufsiot  uia. »  o.stsg  inches 

PIN  FLIGHT’  0.12S9  INCHES 

FIN  THICrHF SS*  0.020?  INCHES 

FIN  SPACING,  CEN Tfc R-TO-CCNT  L.3*  3.0800  INCHES 


tuae  arrangement 

TRANSVERSE  PITCH’  1.0223  INCHES 

LONGITUDINAL  PITCH’  0.83/4  INCHES 

BANK  HEIGHT’  236.0H35  INCHES 

BANK  wlOTH*  136.0704  INCHES 

THETA*  0.65SA 

DELPHI-  -0.3720 


INTEGER  PARAMETERS 

TYPE  CF  CFCSS  FLO  *  ARRANGEMENT*  7 

TYPE  CF  FH  P  <  JP  I  Lc  ’  1 

NUMBER  JF  H G n S ’  4 

NUMBER  OF  PASSES*  2 

VERTICAL  ROWS  CONTAIN  229  TUBES 


PISCELL AVEOUS  VARIABLES 

given  hem  transfer  fate*  iooooojo.qooo  btu/hr 

THERMAL  CCMuJlT IV  ITT  CH  TUJti  MATERIAL*  22  0.0000  3TUMR-FT-F 

THfcRRAc  CONDUCT  IVlTY  OF  c  IN  V.ATERUL-  113.0000  RTU/HR-FT-F 

GAS  CONSTANT*  53.34‘Ol  F  T-L  3F/L  BM-R 


PENALTY  PARAMETERS 

12.C00000000  0.0  0.0 
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CTHER  VALUES 

HEAT  TRANSFER  RATtS 

Q  DO  T  l  *  99SJ2  30.0000  B1U/HR 

QOUT2*  99  9999-1.  jOOJ 

tJL’OTJ*  9<5<;999-*.000,> 

HEAT  TRANSFER  DIFFERENCES 
OELCl*  672J.OJO)  6TU/h< 

DELU2*  -o 713.OOO0 

DELCj*  4.0JUO 

PENALTY*  11.9919 


ccnstpaints 

DIAMETER  RATIO* 

TUBE  THICKNESS* 

TOUCH.'!*  -0.1369 

TCUCHl*  -0.002J 

PRO  FH*  -7  ).  0000  E 

PROFC*  -30.0000  F 


1.4513 

0.013  4  INCHES 


VALUES  TO  BE  OPTIMIZED 
VOLUME*  71.3214  CJ.FT. 

AREA*  10105.4219  iJ.FT. 

A  IR  HP*  92.  7296  HP 


AIR  PRESSURE  DROP* 

tubes ioe  pressure  drop* 


0.0710  PSI 

0.1334  PSI 


OBJECTIVE  FUNCTION*  71.8214 

l  PROGRAM  CALLS  TO  ANA  LIZ 


ICALC  CALLS 

1  l 

2  1150 

3  1 
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APPENDIX  D 

HEDSUP  PROGRAM  LISTING 
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